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ABSTRACT 


In order to obtain a better understanding of the results 
of ohe Mariner IV occultation experiment, a theoretical study 
of tne structure of the neutral atmosphere and ionosphere on 
Inrs has been made. The distribution of temperature in the 
lower atmospnere is obtained from a steady state solution of 
the equation of radiative transfer, allowing for vibrational 
relaxation of the radiating gas, and for the absorption of 
solar energy in the atmosphere. Possible relaxation of the 

vibration - rotation near infrared bands of CO in v.-hich 

o 

solar energy is absorbed is considered, and atmospheric 
uransraissicn i unctions an the near and far irsTrared bands of 
CCp calculated using a modified version of the Blsasser 
band model, taking both Lorentu and Doppler broadening of the 
spectral lines into account. The required input data are 
chosen from, the most recent estimates of the surface para- 
meters and atm.ospheric com.position available at this tim.e. 
Possible effects of CC2 photochemistry on the atm.ospheric 
structure are estimated from, a separate calculation of the 
photochemdcal equilibrium.! distribution of CO2 and fts dis- 
sociation products. Bsti.mates of the solar heating produced 
by this photochemical activity ^re obtained simultaneously 
v;ith the distribution of the gases. From consideration of 
the possible influence of turbulent mixing and molecular 
diffusion in altering the distribution of gases from 
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photochemical equilibritm conditions, upper and lov/er limits 
are derived for the neutral particle number density in the 
observed ionosphere. An ionospheric model consistent v;ith 
the lov;er atmospheric structure is proposed. 

The results of this study indicate that the ruTi of 
atmospheric temperatiu*es on I'hrs at the time and place of 
I‘*a.riner IV occultation v;as colder than had been estimated 
fromi previous thermal calculations v;hich v;ere used as a basis 
for an E-layer interpretation of the observed I-'artian ionos- 
phere • At the samiC time, the thermal structures obtained 
here are considerably v;armer than the heuristic miOdels pro- 
posed to explain the observed ionosphere as a lov;-density, 
diffusion-controlled F2-peak, and demonstrate the unaccept- 
ability of this interpretation. For surface temperatures 
ranging from 175®K to 200°K, surface pressures ranging from 
5 mb to S mb, and an atmospheric composition of to lOOfs 
CO 2 , the tropopause level is found betv;een 2 km and 5 km 
above the surface of I-ars. The convective layer is shallov: 
(less than 5 km thick) in all cases considered, and it is 
estimated that at a ground temperature below ^ 140'^a, no 
troposp!!bre would develop. For surface temperature and pres- 
sure and composition v^ithin the above range, atmospheric 
tem.peratures decrease with altitude to values of 140®K - 150°K, 
13 5 ‘’X - 145 “K, and 130*K - 140“ X at 50 km, 100 km and 125 km, 
respectively. The influence of the absorption of solar near 
infrared and ultraviolet energy on the atmospheric structure 
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belov; ionospheric levels is found to be srnll. Cn the basis 
of the calculated thermal models for the lov;er atmosphere, 
it is concluded that the neutral particle number density in 
the observed ionosphere was betv;een 10^^ mol/cm^ and 
10^^ mol/cm^, about 10^^ mol/cm^ if mixing extends to near- 
ionospheric levels, and about 10^^ mol/cm^ if the atmospheric 
gases are gravitationally separated above the region of maxi- 
mum photochemical activity. Thus, the Lhrtian ionosphere 
observed by I>^riner IV is of the Fl-type. 

It is demonstrated that an Fl-peak of ions in an 
atomic oxygen environment is unlikely, since the requirements 
for suppression of an F2-peak at higher levels can be met only 
by contradiction of a basic assur.TOtion, A peak composed of 
booh CO 2 and 0^ ions in a CO^ environment appears to be a 
more acceptable alternative, and suppression of an F2-peak 
can readily occur in this case if the 0 concentration at the 
Fl-peak is less than one-tenth of the CO^ concentration. It 
is thus suggested that the I'hrtian atmosphere is raxed up to 
near-ionospheric levels, although at and above the observed 
ion maximum gravitational separation may be present. The 
neutral particle density at the observed peak is /v 2 x 10^^ 
mol/cm , predominately C 02 * ^ ^^2 “ ^2 distribution 

consistent v/ith the observed ion profile is calculated from 
basic ionospheric theory. The present ionospheric model 
allows for the development of a true thermospheric structure 
with a positive temperature gradient in the upper L'artian 
atmosphere. 
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INTRODUCTION 


Formulation of the Problem 

The principal problem treated in this paper is the de- 
termination of a model vertical structure of the Martian 
atmosphere which is both physically sound and consistent 
with the most recent observational findings. The basic ver- 
tical structure is obtained from solution of the equation of 
radiative transfer in the presence of energy sources and 
sinks, allowing for vibrational relaxation of the infrared- 
active constituent. 


Motivation and ob.iective . 

The brilliant success of the Mariner IV occultation 
experiment resulted in significant new estimates of the 
s\irface parameters and vertical structure of the lower Martian 
atmosphere, as well as the first observational evidence of 
the existence of a Martian ionosphere. The occultation data 
clearly indicate that the atmospheric pressure, density and 
scale height near the surface of Mars are all considerably 
lower than had previously been supposed. They also indicate 
that the ionosphere is located nearer the surface and has a 


far smaller electron density than had been predicted on the 


basis of classical ionospheric theory. Although free from 
the terrestrial interference which has limited the precision 
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of Earth-based observations of Mars, the Mariner IV data have, 
nevertheless, not fully resolved the problem of the Martian 
atmosphere, and have left room for a variety of interpretations. 
However, the basic indications of the I-lariner IV occultation 
experiment are indisputable and must serve as a guidepost for 
further theoretical and experimental investigations of the 
Martian atmosphere. 

Principal discussion in the literature of the Mariner IV 
findings has been with regard to the identification of the 
type of the Martian ionosphere observed in the occultation 
phase-shift data. The relative importance of the various 
processes responsible for the formation of an ionosphere depends 
principally on the number of neutral particles present locally. 
Thus the type of ionosphere observed by Mariner IV depends 
upon the molecular number density at ionospheric height 
( 'v 125 km). The molecular number density at this level cannot 
be derived from the occultation data and must be estimated 
from theoretical modelling. Hence, a knowledge of the vertical 
structure of the lower Martian atmosphere is essential to an 
understanding of the Martian ionosphere. 

Meaningful models of the vertical structure of a planetary 
atmosphere can be obtained from consideration of basic physical 
principles when certain characteristics of the planet and its 
atmosphere are known. These characteristics include the sur 
face pressure and temperature, the type and amount of infrared- 
active constituents present, and the spectral distribution of 
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solar energy incident at the top of the atmosphere. By speci- 
fying these quantites and by considering in detail the internal 
sources and sinks of energy together with the mechanisms of 
vertical energy transport and the important photochemical 
reactions, an equilibrium height-distribution of temperature, 
density and chemical composition can be determined for the 
planetary atmosphere. Then, from this vertical struct\ire, a 
better understanding of the ionosphere of the planet can be 
achieved. 

Previous theoretical studies of the atmospheric structure 
on Mars are of limited practical value for two basic reasons. 
Firstly, they have dealt with simplified models and neglected 
physical processes which could be important such as solar 
heating, photochemistry, and non-local thermodynamic equilib- 
rium effects. Secondly, they have been based upon assumptions 
of s’orface pressure and atmospheric composition which are no 
longer realistic in view of the most recent estimates of these 
parameters. 

The primary objective of the present study is the deter- 
mination of an equilibrium distribution of temperature, density 
and chemical composition for the Martian atmosphere in accord 
with the latest observations, considering all processes 
important in the energy balance of the atmosphere. 
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Scope of present study . 

From our understanding of the basic physical processes 
which govern the vertical temperature distribution in the 
Earth *s atmosphere, some insight into the energetics of 
the Martian atmosphere may be gained. 

Near the surface of Mars, convective energy transport 
will be large, due to the discontinuous nature of the lower 
boundary, and a troposphere will be established, extending 
upward to the level at which the lapse rate of temperature 
becomes stable against convection. Since surface temperatures 
on Kars are generally lower than are found on Earth, and the 
primary constituent of the Martian atmosphere, C02» is a 
strongly radiating gas, one would expect the Martian convec- 
tive layer to be thinner than the Earth’s troposphere. 

Above the convective layer, radiative transfer by the 
large amounts of CO2 present on Mars will largely control the 
temperature structure. At levels in a planetary atmosphere 
where temperatures are moderate and where there is an abun- 
dance of polyatomic molecules, radiative fluxes will determine 
the distribution of temperature. The radiative transfer 
process in the Earth’s atmosphere is complicated by the 
presence of three infrared-active gases {CO21 0^ and H2O) with 
overlapping absorption spectra. Furthermore, a sizable tem- 
perature maximum is produced in the middle atmosphere of the 
Earth by the large amount of solar energy absorbed in the 
dissociation of O2 and photochemically produced O3. On Mars, 
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however, the problem of radiative transfer is considerably sim- 
plified, since observations indicate that the Martian atmosphere 
is composed primarily of COg* In the presence of large amounts 
of CO2 and moderate pressures ( > 10 ”^ mb), the absorption of 
solar energy would probably have little effect on the local 
temperature, due to the very large radiative fluxes of planetary 
origin* One expects, therefore, no temperature maximum of any 
consequence in the middle atmosphere of Mars, with the temper- 
ature structure largely determined by radiative transfer in 
the 15 ^ CO2 band. 

At higher levels in the Martian atmosphere where pressure 
is low ( < 10 mb), the absorption of solar energy may exert 
a stronger effect on the local temperature, even if CO2 is 
present in abundant quantities. The transfer of energy between 
the molecular kinetic mode and the radiative field via the 
vibrationally excited states of GO^ will be inhibited at higher 
levels, where the radiative lifetime of these excited states 
becomes shorter than the collisional deactivation or "relax- 
ation” lifetime. At these levels, with the linkage between 
matter and radiative field broken by ’’relaxation’*, any external 
source supplying energy directly to the kinetic mode might 
considerably alter the temperature structure. In the Earth’s 
atmosphere , vibrational relaxation in the 15 /f CO2 band occurs 
between 75 km and 35 hm. The vibrational relaxation level on 
Mars should be found at a lower altitude, however, due to the 
lower surface pressure and colder atmospheric temperatures 
anticipated. 



The photodissociation of CO 2 and that of photocheraically 
produced ©2 and 0^ vdll result in the deposition of solar 
ultraviolet energy in the region of the Martian atmosphere 
where vibrational relaxation is occuring. Because of the rel- 
ative magnitudes of the lifetimes for the collisional and 
radiative processes involved, this energy will pass into the 
kinetic mode with an efficiency less than but of the order of 1 
In addition, CO 2 absorbs strongly in the lyif - 6^ near infra- 
red bands where the solar flux is strong also. However, most 
of this energy should be scattered at higher levels in the 
Martian atmosphere where these bands also relax, and passage 
of the solar near infrared energy into the kinetic mode should 
be minimized. Still, since the amount of energy in the near 
infrared solar spectrum is much greater than that found in 
the ultraviolet spectral regions, it is possible that the 
absorption of near infrared energy would exert a larger influ- 
ence on the temperat\ire at higher levels than would the ultra- 
violet absorption. 

In higher regions of the Martian atmosphere, turbulent 
mixing must become insignificant and the constituent gases 
must separate diffusively in the gravitational field of the 
planet. In the terrestrial atmosphere, the transition between 
homogeneous and inhomogeneous atmospheric regions is rather 
abrupt, due to the highly stable conditions which result from 
the sharp rise in temperature above the Earth *s mesopause. It 
is Improbable that a region of such high stability is present 
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in the upper atmosphere of Mars, as temperature gradients should 
be relatively small} therefore no level can be set with any 
definition for the turbopause on Mars* 

In the upper Martian atmosphere, as in the Earth’s upper 
atmosphere, where the vertical temperature gradient becomes 
significant as the temperature rises to its thermospheric values, 
the heat conduction process will begin to play an important role 
in the vertical transport of energy* At altitudes where CO 2 
is largely dissociated, radiative loss by forbidden transitions 
of 0 atoms and vibrational-rotational transitions of CO mole- 
cules must be considered* In this region, equilibrium can be 
assumed to exist between the photoionizing solar energy absorbed, 
the divergence of the conductive flux and the loss of radiant 
energy by 0 or CO emission* Here the temperature structure 
must be obtained from a solution of the heat conduction equation 
with radiative sources and sinks* It is unlikely, however, 
that temperature gradients below the level of the l^artian ionos- 
phere would be large enough to make the divergence of the con- 
ductive heat flux comparable in magnitude to the heating and 
cooling in the \5f4 00^ band, even if CO 2 undergoes moderate 
dissociation below the ionospheric level* 

By a proper synthesizing treatment of the various aspects 
of the atmospheric problem outlined above, a totally self- 
consistent solution for an equilibria structure of the Martian 
atmosphere will be obtained numerically by an iterative solution 



• of the energy balance equation. Not only will the resulting 
distributions of temperature, density and composition be 
self-consistent, but, furthermore, they will be consistent 
with the boundary conditions suggested by the most recent 


observations of Mars* 
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Earth~Based Observations of Ifers 

Since the txirn of the century, numerous experiments have 
been performed to determine the surface and atmospheric condi- 
tions on I^rs, Radiometric, photometric, polarimetric and 
spectroscopic methods have all been employed and a sizable 
body of observational data has been amassed. In addition, the 
observations have been clarified a^d supplemented by many 
theoretical works. 

Estimates of the surface temperature on Mars can be readily 
obtained from infrared and microwave radioraetry, and, conse- 
quently, this parameter has been fairly well established for 
several decades. Interpretation of photometric, polarimetric 
and spectroscopic measurements is far more difficult, and no 
clear picture of the surface pressure or atmospheric composition 
on V^ra emerged from any of these techniques until quite 
recently. Knowledqe of the vertical structure of the Martian 
atmosphere cannot bo gained from observations made from the 
surface of the Earth. Since the major infrared-active consti- 
tuent of the atmosphere of Mars is also an important component 
of the terrestrial atmosphere, interpretation of observed 
spectra of Jfeirs to obtain a vertical temperature distribution 
is not practicable. 

liiose Earth-based observations of greatest historical 

importance and relevance to the problem are briefly reviewed 

, $ 

here. 



The first infrared radiometric measurements of the Martian 
svxrface were made by Coblents and Lampland (1923, 1927) 

Lowell Observatory and by Pettit and Nicholson (1924) at Mount 
V/ilson during the oppositions of 1922, 1924 and 1926. Using 
various filtering techniques and lunar calibration, these exper- 
imenters made latitudinal and longitudinal scans of the planetary 
disk to determine the diurnal and seasonal variations of temper- 
ature as a function of latitude. An extension of this work was 
carried on by Lampland in a systematic program of ^!^rtian radi- 
ometry at Lowell Observatory between 1926 and 1943* A summary of 
all measurements of surface temperatures on 1'fe.rs prior to 1952 
has been given by de Vaucouleurs (1954), along with a detailed 
discussion of experimental methods, the interpretation of meas- 
urements and the problems encountered in making observations of 
this type, 

Hess (1950) constructed an isotherm map for Phrs during 
northern hemisphere winter from measurements made along the noon 
meridian by Coblentz and Lampland (1927). Gifford (1956) ana- 
lyzed the much larger body of data obtained over the years by 
Lampland, and determined the average diurnal temperature vari- 
ation near the Martian equator, the average temperature vari- 
along the noon meridian for the four seasons, and seasonal noon 
meridian isotherm maps for southern hemisphere summer, winter 
and fall. His southern hemisphere summer map agrees in general 
with the map of Hess (1950). 
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.During several weeks following the 1954 opposition, Sinton 
and Strong (I960) obtained drift curve measurements of surface 
temperature on I-2ars at the Mount Paloraar Observatory. Their 

results agree quite well with the earlier observations and are 
summarized in Table I (Rasool, 1963). 

Determinations of the raicrcn-Tave brightness temperature of 
the I4artian disk by Mayer, ^ al. (1953) and Giordmaine, ^ al. 
(1959) resulted in soraev/hat loxirer values of surface temperature 
than those obtained from infrared radiometry. ^kyer (1961) has 
explained this discrepancy as a consequence of the microwave 
radiation originating a few centimeters beneath the surface of 
the planet, where, on the sunlit side, the temperature is low^er 
than that of the surface itself. 

Inking use of the observations of Sinton and Strong (I960) 
and Giordmaine, et (1959) as well as the results of Gifford 
^^956), Mints (1961) has derived the latitudinal variation of 
surface temperature on Mars. Figiure 1, from Mints (1961), shows 
the variation of daily mean, maximum (noontime) and minimum 
(sunrise) temperatures with latitude at northern hemisphere 
v/inter solstice. Although based on a subjective synthesis of 
observational evidence, Figure 1 must depict the variation of 
surface temperature on Mars with reasonable accuracy. 

Atmospheric composition and surface pressure . 

In 1947 at McDonald Observatory, the 1.5’^ and 1.60// 
vibrati^^rotatlon bands of CO were observed in the spectrum 



TABLE I 


TEI'iPERATURE OF THE MARTIAN SURFACE 
(After Rasool, 1963) 


Maximum temperature at equator 

Mean amplitude of diurnal variation 
Noon to sunset 
On Earth in desert 

Night side temperature cannot be measured 
but probably can be estimated at the 
equator 

Day side temperature at poles 

Mean temperature of day side 

Mean temperature of whole planet 


A, 300®K 

Af 60“K 
V 30® K 

/V 200® K 
v220®K 
N 260®K 
N 230®K 



Temperature, In 



Temperoture , in 
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of Kars by Kuiper (1952) (See Figxire 2.). This observation v;as 
confirmed by Kuipcr in the following year when he detected three 
additional CO 2 bands at 1.96y^, 2.01/< and 2.06yy. By comparing 
the intensity of the bands observed in the Kartian spectrum with 
those observed in the spectrum of the sun, he estimated the 
amount of COg present in the atmosphere of Mars to be /v 4.4 
m-atra, about twice the amount found in the Earth’s atmosphere 
(Kuiper, 1952)# 

Grandjean and Goody (1955) redetermined the concentration 
of CO 2 in the i</jartian atmosphere from Kuiper’ s observations, 
taking into account the fine structxare of the band and the shape 
of the rotation lines. They found that for a surface pressure 
of 100 mb, the total amount of CO 2 on Mnrs would be thirteen 
times greater than the amount present in the atmosphere of the 
Earth, or about 30 m-atm. 

In 1963, using high-resolution spectroscopy, Kaplan, Munch 
and Spinrad (1964) at Mount Wilson detected the rotational 
lines of CO 2 in a weak band near 8700 A, and, from laboratory 
determinations of the line strengths, estimated the amount of 
CO 2 in the Martian atmosphere to be 55 1 20 m-atm. 

Early attempts to detect HgO in the atmosphere of Mars 
(Adams and Dunham, 1937; Adams, 194D were unsuccessful and an 
upper limit of 10‘^ g/cm^ precipitable vater was set by the 
investigators (Dunham, 1952). Later negative observations by 
Kiess, et (1957) had resulted in a larger estimate (8 x 10 * 
g/cm^ precipitable water) for the amount of HgO which could be 
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present on I*iars and yet remain undetected. Finally, in 1963, 
several independent observations positively verified the 
presence of H^O in the 14artian atmosphere. Dollfus (1963) 
detected H 2 O bands near 1.4y^ using a polarizing photometric 
method and estimated the amount of precipitable water on Mars 
at 2 X 10"^ g/cm^. Danielson, ^ al. (1964), from their 
spectroscopic measiirements in the 2,7yt( band, suggested that less 
than 4 x 10"^ g/cm^ precipitable water is present, while Kaplan, 
Munch and Spinrad (1964) arrived at an amount of (1.4 - .7) x 10 ? 
g/cm^ from their observation of the rotational lines of H 2 O 
near 6300 A. 

Spectroscopic observations at Mount V/ilson during 1933 and 
1934 by Adams and Dunham, rediscussed by Dunham (1952), failed 
to detect any evidence of O 2 in the ^^artian spectrum, and it was 
concluded, therefore, that the amoxmt of O 2 present could not 
exceed 2.4 m-atm. Richardson (1957) and Kiess, ^ al. (1957) 
confirmed this negative result for 02»and Kaplan, Munch and 
Spinrad (1964), from the absence of O 2 in their spectra, derived 
an upper limit of 70 cm-atm for this gas on Mars, a limit lower 
by a factor of three than that siiggested by Dunham (1952). 

To date, examinations of the spectra of ^^ars for evidence 
of various other gases (Kuiper, 1952; S inton, 1961) have had 
only negative results , and upper limits have been given for the 
amounts of these gases. (See Rasool (1963) for a summary of the 
possible abundances of gases in the I^artian atmosphere.) 
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From measurements of the polarization of sunlight scattered 
by the planet, Lyot (1929) estimated the total pressure of the 
atmosphere on Mars to be less than 25 mb. This estimate was 
superseded by that of Dollfus (1957) who, at the Pic du Midi 
Observatory, made polarimetric measurements across the lartian 
disk at different phase angles and in different colors. By 
several methods of data reduction, Dollfus obtained a value of 
about 85 mb for the surface pressure. De Vaucouleurs ( 1954 ), 
from visual observation of the brightness of Martian markings at 
the Peridier Observatory in 1939, deduced a value of 80 1 13 mb, 
while photometric determinations of atmospheric density 
(Sytinskaya, 1962) have given surface pressures ranging betv/een 
about 60 mb and 110 mb. Kaplan, Munch and Spinrad (1964), by 
combining their estimate of CO2 abundance with spectroscopic 
observations by Kuiper (1963) and Sint on (1963) of the strong 
Martian bands near 2yi(, estimated the total atmospheric presstire 
on Viars to be 25 1 15 mb. This pressure estimate suggested the 
possibility that CO2 is a major constituent of the L'lartian 
atmosphere. 

1 In reviev/ing the various techniques which have been used to 

estimate the surface pressure and atmospheric composition on Mars, 
Chamberlain and Hunten (1965) concluded that both polarimetric 
and photometric methods of pressure determination are unreliable. 
Because both techniques vrould actually measure scattering by any 
submicron particles present in ciCluix^xon to molecular scattering, 
they may be worthless for determining atmospheric pressure. 
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According to Chamberlain and Hunten, spectroscopic observations 
are, in principle, capable of yielding both atmospheric pressure 
and composition. In their view, hovfever, observations made 
prior to 1965 are not refined enough to estimate either para- 
meter v/ith great accuracy. They underlined the need for an 
improved determination of CO 2 abundance from a weak band, 

Belton and Hunten (1966) observed the weak 1.05yii band of 
CO 2 in the spectrum of ^kr3 at Kitt Peak National Observatory in 
1965 and derived an atmospheric abundance of 6S 1 26 m-atm for 
C02> corresponding to a partial pressure of 5 1 2 mb. They have 
coupled this result v/ith Gray’s (1966) analysis of Kuiper’s 
(1963) data for the strong 2^ bands to obtain a sxirface pressure 
in the range of 5 mb to 13 mb, with a preferred value of 6 mb. 
Thus, their results indicate that CO 2 is the dominant consti- 
tuent of the Martian atmosphere. 


/ 


/ 



Pre - Mariner IV Theoretical Investipiations 
of Atmospheric Structure 
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Lower atmosphere , 

A number of theoretical models for the lower ^fe,rtian atmos- 
phere have been constructed from radiative transfer consider- 
ations. Most of the models, however, are physically oversimpli- 
fied and based upon earlier estimates of surface pressure and 
atmospheric composition which are now known to be unrealistic. 

Goody (1957) computed the vertical distribution of tempera- 
ture in the I'lartian atmosphere taking both convective and 
radiative energy transport into account. Assuming a ground 
temperature of 270®K suggested for noontime equatorial regions 
by Hess (1950), a surface pressure of G5 mb suggested by the 
polarization measurements of Dollfus (1957) and a tropospheric 
lapse rate of 3*7®K/km suggested by the small value of the 
acceleration due to gravity on Ifers, Goody obtained the temper- 
ature structure for two different compositions. For one model 
he chose CO 2 as the only infrared-active constituent, with a 
concentration suggested by the observations of Kuiper (1952) and 
the analysis of Grand jean and Goody (1955). For the other model 
he chose H 2 O as the only infrared absorber with a total amount 
of 10 cm precipitable water. In the case of CO 2 alone, the 
tropopause was found near G, 5 km at a temperature of -v 240° K, 
with temperatures falling off with height to an approximately 
constant value of 134®K above 90 km. Since the temperature 
gradient in the region above the tropopause is only just stable 
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in this model, the exact trepepause height is somewhat indef- 
inite. In the case of HgO alone, a 25 km tropopause height was 
obtained, with a tropopause temperature of 17^® K. At higher 
levels in this model the atmosphere is more nearly isothermal, 
with temperatures approaching a constant value of 153* K. Goody 
also fo\ind that decreasing the HgO concentration in this model 
would cause the stratosphere to warm, but that the maximum value 
of tropopause temperature that could bo achieved v/as 227® K, 
with the tropopause at 11.6 km, 

Arking (1962) examined the effects of three assumptions 
commonly used in treating the problem of vertical heat transfer 
in planetary atmospheres* These assumptions are (1) that the 
atmosphere is infinite in extent, (2) that absorption and 
emission are grey and (3) that no convection is present. He 
assumed that the atmosphere is completely transparent to solar 
radiation and that conditions of local thermodynamic equilibrium 
exist at all levels. He found that when the atmosphere is 
optically thick (optical depth 3 or 4) , the temperatures remain 
v/ithin 1^ or 2^ of the results for an infinite atmosphere (the 
Milne problem), but that for optically thin atmospheres (optical 
depth A/ 1 or less) the temperature deviations may be greater 
than IC^. He also found that the effect of a spectral window on 
the cooling of the planetary surface cannot be offset by a band 
of equal strength elsewhere in the spectrum; raising the 
absorption coefficient somewhere in the spectrum has a much 
smaller effect than lowering it, and this non-grey cooling 
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effect becomes much more important the thicker the atmosphere 
is. In calculating a grey-convective model designed to approx- 
imate conditions in the ^fertian atmosphere, Arking assumed an 
effective planetary temperature of 217“K, a total optical thick- 
ness of 0,5 and the decrease of opacity with height v/ith an 
»* absorption scale height” of 17 km. He found a tropopause at 
the 8 km level, a surface temperature of 235®K and a temperature 
of 'v 188® K at an altitude of 30 km, 

Ohring (1963) computed the vertical run of atmospheric 
temperatures on Kars assuming a surface temperature of 230® K 
and a surface pressure of 85 mb, with a tropospheric lapse rate 
of 3*7®K/km and a CO2 concentration of 2^ by volume. He obtained 
a tropopause height and temperature of 9 km and 196® K, respec- 
tively. Above the tropopause the temperature decreases at a 
rate only slightly less than adiabatic to the rather lov/ value 
of about 90® K near 42 km, Ohring comments that at this temper- 
ature CO2 would condense and that this condensation phenomenon 
might have some connection with the blue haze observed in the 
Kartian atmosphere. 

More meaningful models of the atmospheric structure on Infers 
were developed by Prabhakara and Hogan (1965), who were motivated 
by the fact that none of the earlier models described above had 
taken the absorption of solar energy in the atmosphere into 
account, as well as the fact that the input values of the para- 
meters employed in these models had become outdated. Prabhakara 
and Hogan reexamined the radiative equilibrium temperature 
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structure of the I4artian atmosphere, including the absorption 
of solar ultraviolet and visible energy by oxygen and ozone 
and solar near infrared energy by CO^ in the 1^ - 6/f bands 
in the calculation of atmospheric heating. They theoretically 
calculated the transmission functions of CO , making use of 
the strong line approximation to the ♦’statistical” model to 
simulate band absorption, and used these transmission func- 
tions to evaluate the absorption of solar energy in the near 
infrared and to investigate the radiative transfer in the far 
infrared. For various spectral intervals in the near and far 
infrared, they derived the theoretical band parameters involv- 
ing the line intensity and mean ratio of line half -width to 
line spacing from the transmittance tables of CO 2 presented 
by Stull, (1963). 

Prabhakara and Hogan used the basic photochemical theory 
of 0^ production to determine a vertical distribution of 0^ 
consistent with the computed radiative equilibrium temperatxare 
structure. They integrated the equation of radiative transfer 
numerically to obtain the temperature distribution, avoiding use 
of the empirical relationships commonly involved in the pressure 
dependence of CO 2 absorption, and also calculated the infrared 
flux transmittance without any simplifying assumptions. 

Instead of a ”time marching” approach to the radiative transfer 
problem, in which a final solution requires the rate of heating 
to become vanishingly small, they used a steady state approach^ 
requiring equality between absorbed and emitted energies for 
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all levels at each step. They considered both radiative and 
convective energy transport in their solution of the radiative 
transfer equation; whenever the calculated lapse rate tended 
to be superadiabatic it was replaced by the adiabatic lapse 
rate. 

In their study, Prabhakara and Hogan calculated radiative 
equilibrium temperatures from the surface of TAars to the 100 km 
level for six models with surface temperatures ranging from 
230 ® K to 270®K, surface pressures from 10 mb to 50 mb and CO^ 
concentrations from 4 ^ to 44^ by volume. In all of their 
models, they found the tropopause at levels below 10 km, v/ith 
the temperature above the tropopause steadily decreasing to a 
value of about 155** K in the upper layers* Their results indi- 
cated definitely that no temperature maximum is produced in 
the i'-a.rtian atmosphere by the absorption of solar energy in 
the ultraviolet and visible by 0^ and 0^ or in the near infra- 
red by CO 2 . The maximum 0^ number density occurred at the 
surface of V&rs with a gradual decrease with altitude, and the 
total amount of ozone present was found to be about one-tenth 
of the amount found in the Earth*s atmosphere ( 0.3 cm-atm) , 

The total ultraviolet energy absorbed in the Martian atmosphere 
by 0^ and 0^ was found to be comparable to the near infrared 
energy absorbed by However, the vertical distribution 

of absorbed energy showed that below n/ 20 km, Og and 0^ 
aosorption is comparable to GO absorption, while above this 
level CO 2 absorption becomes much larger. 
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Of the six models presented by Prabhakara and Hogan, their 
Model I most closely approximates conditions in the lower atmos- 
phere of Mars as revealed by Mariner IV. This model was widely 
referred to in post-Mar iner IV studies of the I'fertian atmosphere, 
being directly adopted as a lower atmospheric model by some 
(Chamberlain and McElroy, 1966) and conveniently ignored by 
others (Johnson, 1965; FJeldbo, ^ al. , 1966) in constructing 
ionospheric models. In ^fodel I of Prabhakara and Hogan, consid- 
ered an extreme case at the time, a 10 mb surface pressure was 
adopted along with a ground temperature of 230* K and a CO 2 con- 
centration of by volume. The remainder of the atmosphere 
was assumed to be N 2 * The complete vertical structure of the 
atmosphere resulting in this case is shown in Figure 3. The 
tropopause appeared near the 3 - 4 km level at a temperature of 
A' 212*K, with a decrease in temperature above this level to a 
constant value of 156,5*K at the top of the model. The pressxire 
obtained at the 100 km level was ^ 1,2 x 10"^ mb, and a number 
density of ^ 5 x 10 mol/cxar at this level resulted. 

Another result of the study of Prabhakara and Hogan indica- 
ted that solar heating in the near infrared bands of CO 2 might 
v/arm the Martian atmosphere at some levels by v 20* K above its 
radiative - convective equilibrium temperat\ire structure. This 
effect is illustrated in Figure 4 (an amended version from 
Prabhakara and Hogan (1965)), for a case of 30 mb surface 
pressure and 9j5 by volume CO 2 concentration. In this particular 
case, excluding the absorption of solar energy lowers the temper* 
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Figure 4. Equilibrium temperature profiles calculated for I-Iars:^ 

radiative, radiative-convective and radiative-convective 
with solar heating (After Prabhakara and Hogan, 1965). 
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ature considerably, and the cumulative effect of the lower 
temperatures causes a significant reduction in number density 
at the 100 km level* In the work of Prabhakara and Hogan (1965) 
the existence of local thermodynamic equilibrium conditions 
throughout the atmosphere was implicitly assumed, both in the 
near and far infrared bands of C02* Also, in their calcu- 
lations they did not give consideration to the Doppler broaden- 
ing of the spectral lines or photodissociation of CO 2 . V/ell 
aware of these limitations , they set a limit to the region of 
validity of their results. In the case of Model I, their 
results are admittedly of questionable accuracy above 65 km 
altitude. 

Upper atmosphere . 

No observational evidence of the characteristics of the 
upper atmosphere of I‘'5ars was available prior to the I'ariner IV 
occultation experiment. Upper atmospheric models had been 
proposed, however, from Earth-analogy and energy balance 
considerations involving the solution of the radiative transfer 
and heat conduction equations. 

Adopting a lower atmospheric model from Goody (1957), 
Chamberlain (1962) derived a temperature structure for the 
Martian mesosphere and thermosphere taking into account the 
various heating and cooling and energy transport processes 
involved. He concluded that the CO, present in the upper 
atmosphere as a result of CO^ photodissociation, has a large 
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thermostatic effect, keeping thermospheric temperatures on Mars 
from exceeding /v 1100® K, so that Mars might retain 0, 
Chamberlain obtained a pronounced temperature minimum of 
A/ 76“ K at the Martian mesopause near a/ 130 km altitude, due 
to large radiative cooling by CO^ in this region, and he sug- 
gested a possible link between this low temperature (at which 
CO^ would freeze out) and the I-'feirtian blue haze* He estimated 
that the solar energy absorbed in molecular photodissociation 
is readily radiated away by CO2 and that the thermosphere is 
produced entirely by photoionization energy* This photoion- 
ization energy is lost by a downward flux of heat by conduc- 
tion, coupled with CO radiative losses near from 

excited vibrational states, so that an equilibrium temperature 
distribution results* The thermospheric model obtained by 
Chamberlain is shown in Figure 5, from Rasool (1963)* (In 
this figure, the temperature structure in the lower atmosphere 
(insert) is basically that computed by Arking (1962)*) 
Chamberlain also speculated on the structure of the 
^^artian ionosphere, estimating smaller ionization densities 
than are found in the Earth* s ionosphere* He suggested an E 
region consisting of two distinct portions, X-rays forming 
an E2 maximum at 200 km, and ultraviolet radiation forming an 
El maximum near the mesopause at 130 km* He also anticipated 
an Fl ionospheric region near 320 km, with a noontime electron 
number density of ^ 10^ electrons/cm^ * He concluded that an 
F2 region on ^^ars would not develop to the extent that it does 
on Earth and all ionospheric regions ndght disappear during the 


night* 



Exosphere 
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Figure 5 , Thermospheric model of Chamberlain (1962) and 
Arking (1962) (After Rasool, 1963), 
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McElroy, et al » (1965) also exananed the structure of the 
upper 1‘fertian atmosphere, constructing a number of models vdth 
various constants and boundary conditions in the heat conduc- 
tion equation. They foxmd that CO vibrational cooling is not 
as important as had been suggested by Chamberlain, but that 
the rotational levels of CO and the fine structure levels of 
the ground state of 0 may cause appreciable radiative loss. 

They also found that changing the composition at the meso- 
pause level only moderately affects the temperatxire structure 
as does increasing the total number density at the mesopause 
level. Exospheric temperatures ranged between <v 725® K and 'V 
950® K in the three models illustrated, McElroy, ^ al, also 
found that it appears to be a good approximation to neglect 
radiative losses to space at optical depths greater than unity. 
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Mariner IV O ccultation Experiment 

The highly successful mission of Mariner IV provided inter- 
esting new information on the physical properties of the planet 
Mars. On November 2 B, 1964, the Mariner IV spacecraft was 
launched into orbit about the sun on a trajectory which would 
result in rendezvous with Mars and occultation of the spacecraft 
by Mars as viewed from Earth. On July I5, I965, with data 
acquiring systems operating at their maximum capability, 

Mariner IV approached to within 6IIS miles of the Martian surface, 
passed behind the planet and reappeared again on the opposite 
limb. The encounter sequence and spacecraft trajectory near 
encounter are depicted in Figure 6, after Anderson (I965). At 
the time of the encounter Mars was 2,l6 x 10^ km from the Earth 
and 2.32 X 10 km from the sun with an Earth - sun - Mars angle 

of 64®. The subpolar point at encounter was 15.3® N of the 
Martian equator. 

Magnetometer and high and low energy particle detector 
measurements in the vicinity of Mars all indicated that only a 
very weak magnetic field { < 10”^ gauss) can be present (Smith, 
et al., 1965; O’Gallagher and Simpson, 1965; Van Allen, et ^. , 
1965), and photography of the planet revealed that the Martian 
surface, like that of the moon, is heavily cratered (Leighton, 

£t al. , 1965). Of particular relevance to the present study 
is the Mariner IV occultation experiment (Kliore,^ ^. , 1965) 
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I'ariner IV encounter sequence and trajectory 
(After Anderson, 1965 )« 



33 


which resulted in estinates of the gross properties of the 

lower atmosphere of Mars, and clearly indicated the presence 
of an ionosphere* 

Experimental techniques . 

The possibility of using the radio waves propagated 
between the Earth and a spacecraft whose trajectory passes 
behind a planet to obtain certain properties of the planetary 
atmosphere has been discussed by Fjeldbo and Eshleraan (1965) 
and Fjeldbo, ^ (1965). As the spacecraft approaches the 

limb of a planet, the phase path, frequency and amplitude of 
the telemetry signal will be altered by the presence of a 
planetary atmosphere and/or ionosphere. The local index of 
refraction in the neutral atmosphere and ionosphere varies 
from unity, so that in these regions, the phase and group 
velocities will depart from their free space values. In 
addition, the mean refractive index along the ray path varies 
with the distance of the path from the planetary limb, and 
bending of the rays from a straight line will occur, along 
with focusing and defocusing of the signal. The changes in 
the phase path, frequency and amplitude of the radio signal 
are thus a function of the distribution of the index of 
refraction in the planetary atmosphere, and if these changes 
are carefully monitored as the spacecraft passes behind the 

planet, the distribution of the index of refraction can be 
inferred. 
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The observed alteration in the characteristics of the 
signal will in general be due to the combined effects of the 
neutral atmosphere and ionosphere of the planet. Since the 
contribution of the neutral atmosphere to the phase path 
change is directly proportional to the integrated atmospheric 
density along the ray path, this contribution will be largest 
when the spacecraft is close to the limb of the planet. On 
the other hand, the effects of the ionosphere on the phase 
path change are proportional to the integrated ion content 
along the ray path, so that these effects will predominate 
when the spacecraft lies at a distance from the limb approx- 
imately equal to the height of the ionosphere above the 
surface of the planet. Since it is reasonable to expect some 
separation in height between the neutral and ionized regions 
causing significant changes in the phase path of the signal, 
it should be possible to obtain information about both 
atmospheric regions, using a single radio frequency. 

The phase changes due to the presence of the planetary 
atmosphere can be obtained by subracting all predictable 
phase shifts due to light transit time, spacecraft motion. 
Earth rotation, etc., from the total phase change during the 
period of the experiment. Extreme accuracy is necessary here; 
in the case of the Mariner IV occultation experiment, the 
total phase change due to all causes other than the planetary 
atmosphere had to be known to 1 part in 10 ^^ (Kliore, et al. , 
1965 ) . Of prime importance is an exact knowledge of the 
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spacecraft trajectory and range rate. If this accuracy is 
achieved, the residual phase changes can be used to estimate 
the properties of the atmosphere and ionosphere by model 
fitting, or alternately, by integral inversion. 

It is obvious, that no unique spatial distribution of 
refractive index corresponds to a given distribution of 
residual phase change as a function of distance from the 
planetary limb. Reason, therefore, must govern the choice of 
an atmospheric model. Fjeldbo, et al. (1965) point out that 
an assumption of a spherically symmetric atmosphere permits 
use of the numerical integral inversion technique in the 
derivation of the refractivity distribution. As is well known, 
however, an ionosphere is not spherically symmetric, with max- 
imum electron-ion number densities being found near the subsolar 
point. A test of the validity of the symmetric ionosphere 
assiuTiption was made by Fjeldbo (I964) and a similar test was 
performed in the present study. Both examinations indicate 
that the non-spherical nature of an ionosphere affects the 
integrated refractive index along a ray path only slightly, and 
that at high latitudes in winter, where the solar zenith angle 
is quite large and varies only over a small range during the 
day, the effects of non-sphericity are negligible. 

I'jeldbo, et al. (I965) also found that changes in phase 
path due to the variation of the electron density in the Earth’s 
atmosphere during the short time period of the experiment should 
be of no consequence and can be excluded from consideration. 
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Results and interpretation of experimenters * 

The estimates of the atmospheric properties on l^Iars 
which resulted from the Mariner IV occultation experiment 
were surprisingly unlike those which had been given on the 
basis of earlier theoretical and experimental work. The 
Kartian atmosphere was found to be colder, more rarified 
and composed of heavier molecules than had been predicted, 
and the ionosphere was discovered to be under -dev eloped 
and located at a lower altitude than had been anticipated. 
These basic atmospheric features were well established by 
Mariner IV. Various atmospheric models have since been 
proposed, however, to clarify the I'lariner observations and 
to give further understanding of the Martian atmosphere and 
ionosphere. These are treated in the next section. 

At Mariner IV entry into occultation, the point of 
tangency of the Earth-spacecraft line on the surface of Mars 
was located at 55® S latitude and 177® E longitude, while 
at exit, the point of tangency was at 60“ N latitude and 
34“ \'I longitude. At the entry point, local time was afternoon 
(13:30 h) in winter, and the solar zenith angle vras a^ 70“; 
at the exit point, local time was near midnight in summer, 
and the solar zenith angle was /v 106“. During entry into 
and exit from occultation, phase, frequency and amplitude 
measurements of the S band (2300 Mcy) Mariner IV telemetry 
signal were made at four receiving stations, two (Pioneer 
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and Echo) at Ooldstone, California and two in Australia at 
Tibinbilla and Woomera. Entry data were obtained with the 
frequency reference of the Mariner IV transmitter provided 
by a frequency standard on Earth, while some of the exit 
data were taken with a crystal oscillator on board the space, 
craft providing the frequency reference. The phase measure- 
ments made during exit using the oscillator as a frequency 
reference lacked precision, and, consequently, emphasis was 
placed on the data obtained during entry into occultation by 
Kliore, et al. (1965). The effects due to the neutral 
atmosphere and ionosphere can clearly be distinguished in 
Figure 7 from Kliore. « al. (1965). In this figure, which 
represents the total phase change on entry into occultation 
as a function of time (from the Goldstone-Pioneer data), the 
maximum ionospheric effect of 10 cycles appears at about 
2:30:10 UT and the final upswing caused by the neutral atmos- 
phere begins at about 2:30:50 UT and reaches a maximum of 

30 cycles at signal extinction (2:31:11,2 UT). 

By relating the maximum phase change {29t 2 cy) , maximum 
frequency change (5.5 - 0.5 cy/sec) and refractive gain at 
the time of signal extinction (1.5 db to 2.0 db) to the 
suriace refractivity and scale height of a model atmosphere 
where density decreases exponentially with height, Kliore, 

H al. (1965) obtained values of 3.6 1 0,2 N units for the 
surface refractivity and 9 1 1 km for the scale height in 



PHASE DIFFERENCE (RESIDUAL SUM) 
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Relation of total phase change to time. Abscissa is UT 
15 July (After Kliore, ^ , 19o5). 
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the lower 30 km of the Martian atmosphere near the entry 
point of tangency. Kliore, ^ (I 965 ) suggest that the 

final ray paths may have grazed an elevated surface feature, 
so that the mean surface ref activity may be slightly larger 
than the above value* 

The inference of density and pressure from the refrac— 
tivity and scale height estimates involves the assumption 

atmospheric composition, since the density ~ refract ivity 
relationship depends upon the type of molecules making up the 
gas. Spectroscopic evidence, coupled with the low scale height 
measured by Mariner IV, indicates that the Martian atmosphere 
consists primarily of CO 2 . Kliore, ^ al. (1965) obtained 
estimates of mass density corresponding to a surface refractiv- 
ity of 3*6 « 0.2 N units for three atmospheric models differing 
somewhat in composition but all having CO 2 as a major constit- 
uent. Temperatures near the surface of Mars were derived from 
the assumed composition (mean molecular weight) and the 
measured scale height of 9 - 1 km. From the values of mass 
density and temperature near the ^^artian surface so derived, 
estimates of surface pressxxre were obtained using the equation 
of state. The surface pressure estimates of Kliore, et al. 
(1965) lie within the range of 4.0 mb to 7.0 mb, in good agree- 
ment v/ith the values obtained spectroscopically by Belton and 
Kunten (1966). Table II from Kliore, ^ al. (1965) contains 
the values for surface temperature, number and mass density 





TABLE II 

SUMMARY OF OCCULTA TICN EXPERIMENT 
(After Kllore, 1965) 


Surface refractivity 


Atmosphere 

3,6 * 0*2 N units 


Scale height 


6 to 10 km 


Sxxrface number density 
100;$ COo 
Up to 20?i 
50^ A 


r/o 


A or N 2 > ^ mixture 


1.9 

2.1 

2,5 


t 0,1 X 10 mol/cm3 
t 0.2 X lOV mol/cm3 
t 0,15 X 10l7 mol/cm3 


Surface mass density 

loo;, CO2 

Up to 20^ A or No. or a mixtxire 
50:5 A 


1,43 * 0.1 X 10*5 g/cm3 
1,5 ± 0.15 X 10-5 g/cm3 
1,75 i 0.10 X 10-5 g/cm3 


Temoerature 
ibo;^ CO2 
Up to 2055 
50 ^ A 


A or N 2 , 


or a mixture 


ISO 1 20’’ K 

175 1 25 *K 

170 t 20 * K 


Surface pressure 
100^ COo 

Up to 20^ A or No, or a mixture 
5 O 5 S A 


4.1 to 5.7 mb 
4.1 to 6.2 mb 
5.0 to 7.0 mb 


Ionosphere 

tjaxiraum electron density (0 *■ 70®) 9 * 1.0 x 10^ el/cm3 

Altitude of maximum 120 to 125 km 

Electron scale height above maximum 20 to 25 km 
Temperature < 200® K at 120 to 200 km 
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and pressure for each of the three models considered by the 
experimenters • 

Measurements made after emersion from occultation 
indicate a surface pressure higher by 50^ than the values 
obtained from the entry data (Fjcldbo, ^ 1966a). 

This discrepancy may be due to the presence of local surface 
features on Mars or due to the departure of the mean surface 
of the planet from a geopotential surface. The temperature 
and lapse rate inferred from the exit data are also higher 
(Fjeldbo, ^ al« , 1966a) as might be expected since the 
emersion point is located in the summer hemisphere. 

Inversion of the ionospheric phase path data obtained 
during entry into occultation indicate a peak electron number 
density of .v 9 x 10^ electrons per cm^ at about 125 km above 
the ^^rtian surface, with an electron scale height above this 
peak of 20 - 25 km. Kliore, et al « (1965) concluded from the 
low altitude of the electron peak and the small electron 
scale height above the peak, that the I-Iartian atmosphere is 
considerably cooler than had been previously anticipated. 

From the absence of a detectable ionospheric phase change in 
the exit data, Kliore, et (1965) set an upper limit on 
the night ime ion density in the region probed by the signal 
at 20 times less than the daytime number density observed, 
or /V 4.5 X 10^ electrons/cm^. A minor ion peak near 95 km 
altitude was suggested by the entry data, which ^x5C ShcVT 



evidence of asynunetry in the ion distribution below the 
peak at 125 km. Kliore, et (1965) have stated that 
their estimate cf electron number density at the peak is 
correct to within lO^^. The ionospheric data derived from 
the occultation experiment are presented in Table II. A 
profile of electron number density vs. height deduced from 
the occultation data by Fjeldbo, al. (1966a) is 
illustrated in Figure d. 
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Post - Mariner IV Atmospheric Models 

The nature of the ionospheric layer detected by Mariner IV 
is the principal problem which has emerged from the results of 
the occultation experiment, and various interpretations have 
been forthcoming in explanation. Closely related to the ionos- 
pheric problem is the question of the structure of the atmosphere 
below. Some of the proposed ionospheric interpretations place 
severe restrictions on lower atmospheric conditions and neces- 
sitate the adoption of atmospheric models which appear to be 
unreasonable on physical grounds. Other interpretations have 
been based upon lower atmospheric models with considerable 
physical backing, and these would seem to be, in general, more 
acceptable. 

In all of the proposed explanations of the Martian ionos- 
phere, analogy with the different regions of the terrestrial 
ionosphere has been widely used. Since the regions of the 
Earth’s ionosphere can be classified in various ways, some basic 
convention must be employed to avoid confusion. Such a conven- 
tion has been suggested by Hunten ( 196 ?) who formally proposed 
categorizing the ionospheric regions according to the physical 
processes governing ion production, redistribution and loss. 

His suggested definitions are: E-laver : A region ionized by 

radiation with an ’’effective” cross-section (i.e, weighted 
by the fractional abundance of the absorber) of a/ 10 cm , 
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where the loss rate is limited by the dissociative recombination 
of molecular ions. Fj^layer: A region ionized by radiation 

with an effective cross-section of /v 10"^7 ^.^2^ 

be limited either by dissociative recombination or by charge 
transfer. F 2-layer : A region of atomic oxygen ions forming a 

peak of electron density above the FI maximum of production rate. 
Here, recombination is limited by the rate of charge transfer 
to molecular ions; as this rate decreases with altitude, dif- 
fusion of ions becomes more important and an ion peak is formed. 
Atomic oxygen is the most abundant constituent in this region. 

As each of the published models may be placed in one of these 
three categories, these general definitions will be adhered to 
in the present discussion* 

F2-laver interpretation . 

It was found necessary by Johnson (1965) and Fjeldbo, ^ 
(1966a), from analogy with the terrestrial ionosphere, to iden- 
tify the ionospheric layer observed by Mariner IV as an F2 ion 
peak. Rishbeth and Barron (I960), in their study of the rela- 
tive importance of recombination and diffusion, demonstrated 
that an F2-type peak occurs very close to the level where the 
loss rate ^ and the diffusion rate (D/H^) are equal, ^ 
where D is the ambipolar diffusion coefficient for electrons 
and ions in atomic oxygen and H is the scale height of atomic 
oxygen, the medium of diffusion. Above the F2 peak, the plasma 
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can be assumed to be in diffusive equilibrium and the plasma 
scale height related to the temperature, if appropriate assump- 
tions are made. By applying the relationships derived by 
Rishbeth and Baron (I960), Johnson (1965) and FJeldbo, et al, 
( 1966 a) concluded that the neutral nximber density at the ionos- 
pheric (F2) peak on Mars must be close to /v 10^ particles/cm • 
Since a number density of -v 10 particles/cra-^ at the surface 
of Mars is indicated by the occultation experiment (Table II), 
a decrease of eight orders of magnitude in particle number 
density within 125 kra is implied by the F2 assumption. Thus, 
the mean scale height below the ionosphere on Mars could be 
no greater than 6,5 km, and the lower Martian atmosphere 
would have to be very cold. Even if, for example, the mean 
molecular weight of the atmosphere were as large as that of 
C0«, the mean temperature would still have to be ^ 130”K, 

By assuming that the atmosphere at ionospheric levels is 
composed almost entirely of 0 atoms, Johnson (1965) and 
Fjeldbo, et (1966a) obtained atmospheric temperatures of 

B5*K and 80®K respectively at ionospheric height from the 
plasma scale height indicated by Mariner IV (Table II). 

Since the plasma scale height was announced by Kliore, et 
( 1965 ) to be fairly constant up to a/ 200 km, Johnson (1965) 
and Fjeldbo, ^ (1966a) concluded that the 60®K - 65® K 

temperature also applied up to the 200 km altitude. They 
found, therefore, that the upper atmosphere of Mars is very 
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cold and nearly isothermal. 

The two thermal structures obtained may be described as 
follows: Johnson (1965) assumed a surface temperature of 

210® K with a convective layer 14 km thick adjacent to the 
ground. The tropopause temperature in his model is 140®K; 
above the tropopause the temperature decreases along the 
vapor pressure curve for dry ice to a value of 85® K at 100 km. 
Above this altitude the atmosphere is isothermal, Fjeldbo, 
et al, (1966a) assumed a surface temperature of 180®K with 
temperature falling gradually with height to the low value of 
50® K near 100 km. Higher up, the temperature rises again, 
and the atmosphere becomes isothermal at ^ 80® K near 135 km. 

There are numerous objections to the two models described 
above, the most serious of which is with regard to the extreme- 
ly cold atmospheric temperatures they require. It has been 
shown convincingly by Prabhakara and Hogan (1965) that radia- 
tive transfer in the I^artian atmosphere is sufficiently strong 
to maintain atmospheric temperatures well above the values 
proposed by Johnson (1965) and Fjeldbo, et al, (1966a), The 
The radiative transfer process has been completely neglected 
in the construction of the above models, which constitute, 
in essence, interpolation between the surface and ionosphere, 
Koreover, careful studies of the Martian thermosphere by 
Chamberlain (1962) and McElroy, et al, (1965) have resulted in 
positive temperatiare gradients in the upper atmosphere in 
every case considered, with exospheric temperatures nearly an 
order of magnitude higher than those proposed by Johnson (1965) 



or Fjeldbo, et al, (1966a). (Hunten (1967) has argued that 
the thermal calculations of Fjeldbo, et (1966b), which 
resulted in an exospheric temperature of /v 130® K, were 
performed using an unstable method, and are meaningless.) 
Therefore, on physical grounds, the F2 interpretation is 
questionable, because of the unrealistically low temperatures 
demanded throughout the atmosphere. 

Gross, et (1966), in proposing an F2 model ionosphere, 
have pointed out that higher temperatures in the Martian upper 
atmosphere are consistent with the >kriner IV observations if 
a positive temperature gradient is present. In the presence 
of such a gradient, the actual scale height can be much 
larger than the isothermal scale height. They obtain an ex- 
ospheric temperature of 550 i 150®K, the large uncertainty due 
to the uncertainties in the heating efficiency and extreme 
ultraviolet solar flux. However, Chamberlain and McElroy (1966) 
have brought attention to an internal inconsistency in the model 
of Gross, ^ a^. (1966), arising from the fact that for the 
atmospheric structure given by this model, the F2 layer would 

be formed well above 125 km. 

In constructing their atmospheric models, it was implic- 
itly assumed by Johnson (1965) » Fjeldbo, et a^. (1966a, 1966b) 
and Gross , ^ ^. ( 1966) that the electron and ion temperatures 
are equal to each other, and that both are equal to the neutral 
gas temperature. However, as pointed out by McElroy (1967) 
and Stewart (1967), it is well established by both observations 
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and theory, that the electron temperature in the Earth*s 
ionosphere exceeds the ion temperature by a factor of 2 or 3 
during the sunlit part of the day, and that at high enough 
levels both exceed the neutral gas temperature. Thus, the 
neutral gas temperatures derived from the observed plasma 
scale height by Johnson (1965), Fjeldbo, ^ al. (1966) and 
Gross, et (1966) are probably all too high. Therefore, 
even harsher conditions would be required by the F2 inter- 
pretation if allowances were made for elevated electron 
temperatiares. 

Finally, a point which would affect the F2 interpretation 
has been raised by Hunten (1967), who contends that double 
*^^^^6i^®r^'tiation of the phase path data to obtain the variation 
of electron scale height with altitude involves accentuation 
of small random and systematic errors. As a result, he finds 
that the observed phase path changes, which were interpreted 
by Kliore, £t al. (1965) as implying a constant electron scale 
height above the electron peak, are actually consistent with 
a scale height gradient as large as 0.6 kn/km. Thus, his 
analysis indicates that rather large scale heights could be 
present some distance above the peak, at an altitude of ^ 200 knu 

These many difficulties with this interpretation appear 
to rule out explanation of the Martian ionosphere as an F2 
region. The fundamental objection to this interpretation is 
that there is no sound physical basis to any of the models 
proposed, and that the low temperatures required to produce 
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an F2 layer at 125 km are impossible to achieve on physical 
grounds. An implication of the xmacceptability of the F2 
ionospheric model is that the derivation of a temperature 
from the ion distribution above the peak is meaningless, 
since the ion distribution must be governed by local produc- 
tion and loss mechanisms, rather than diffusion. 

E-layer interpretation . 

A sounder approach to the problem of the Martian ionosphere 
was made by Chamberlain and McElroy (1966) who combined an 
existing model of the lower atmosphere with their own calcu- 
lations of the thermospheric structure to obtain the neutral 
particle density at the ionospheric level. As a lower atmos- 
pheric structure, Chamberlain and McElroy (1966) adopted 
Model I of Prabhakara and Hogan (1965) (shown in Figure 3) up 
to an altitude of a/ 50 km. They extrapolated this model to 
ionospheric levels, and constructed a model thermosphere upon 
it, treating the various terms involved in the energy balance 
of the higher atmosphere in the manner developed previously 
by Chamberlain (1962) and McElroy, e^ al. (1965)* Strong 
dissociation of CO^ had been assumed in these earlier studies, 
with thermospheric cooling by 0 and CO being important factors 
in establishing the thermal structure of this region. However, 
Chamberlain and McElroy (1966) argue that, in order to maintain 
the generally cool thermosphere indicated by Mariner IV, CO2 
cannot be strongly dissociated, but must exert a strong 
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thermostatic effect, with CO^ radiative losses keeping the 
overall temperature at higher levels low* In addition, they 
argue that due to the long photodissociation lifetime 
( A/ 1 month) of CO^ in the upper Martian atmosphere, mixing 
should control the distribution of gases, so that large quan- 
tities of CO^ must be present at high levels. 

In spite of the assumptions that CO2 is only slightly 
dissociated and that the entire atmosphere is thoroughly mixed, 
both of which have a tendency to produce lower temperatures in 
the upper atmosphere from the radiative cooling of C02» 
Chamberlain and KcElroy (1966) nevertheless find a buildup of 
temperature with height. As a result of the balance betv^een 
solar ionizing the dissociating radiation on the one hand, and 
infrared radiative losses coupled with downward transport of 
heat by conduction on the other, they obtain a true thermosphere 
in v/hich temperatures increase to a value of 400® K or more 
in the exosphere, 

Extrapolating Model I of Prabhakara and Hogan (1965), 

Chamberlain and McElroy (1966) obtain the same number density 

at the 100 km level ( a/ 5 x 10^^ particles/cm3) as ^^)as 

obtained by those authors. At the ionospheric level (125 km), 

3 

they obtain a number density of /v 5 x 10 particles/cm , and 
therefore interpret the observed layer of ionization as a 
Chapman type E-region produced by X-rays at the level where 
the optical depth in the ionizing frequency is unity. Since 
the sun was at a zenith angle of ^ 70® at the time of the 
obsex*vation and since the scale height at the level of the 
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ionosphere must be of the order of 10 km, a neutral particle 
number density of ^ 5 3C 10^ particles/cni^ in the ionosphere 
implies a mean ionization cross-section of the order of 
10“19 cm^* Thus, according to the definitions suggested by 
Hunten (1967) » the ionosphere must be of the E type. For an 
estimated X-ray photon flux of /v 5 x 10 photons/cm sec, and 
a peak electron density of 10^ electrons/cm^ , the effective 
recombination coefficient must be /v 2,5 x 10"^ cm^/sec, a 
reasonable value for an ionospheric region in which dissociative 
recombination limits electron loss, Chamberlain and McElroy 
(1966) argue that although ionization of all constituents occurs, 
various ionic reactions would likely result in most of the 
ionization appearing in the form of NO or 0^ , since these 
species have the lowest ionization potentials of all the con- 
stituents which could be present locally in significant amounts. 

If the I-feirtian ionosphere is an E region, then somehow, 
suppression of FI and F2 regions must be explained, Chamberlain 
and KcElroy (1966) point out that three conditions must be met 
for suppression of an FI region at higher levels: 1) the 

effective recombination coefficient of CO2 must be at least 
10“^ cm^/sec, 2) there cannot be a fast reaction transforming 

CO ■*“ into another ion which has a slower recombination rate, 

2 

and 3) O/CO2 and O2/CO2 abundance ratios must each be less 
than /V 10^, so that 03"^ does not become the dominant ion. To 
satisfy the last condition, there must be little CO2 dissoci- 
ation above 100 km and little diffusive separation at the FI 
level ( ^ 190 km), Chamberlain and McElroy (1966) reason 

C. 
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that v/ith the small O/CO 2 ratio required to explain the absence 
of a detectable FI region, it is unlikely that 0'*’ could become 
the dominant ion below the level at which diffusion becomes 
the main process for ion removal* Thus, no F2 region will 
appear at higher levels* 

Additional thermal calculations for the lower and upper 
atmosphere of Mars have been performed by McElroy (1967) who 
obtains essentially the same atmospheric structure as was 
obtained from the Prabhakara and Hogan (1965) - Chamberlain 
and McElroy (1966) model* A neutral particle number density of 
5 X 10^^ particles/cm^ at the ionospheric level results here 
also, and the E-layer interpretation is adopted. Furthermore, 
the same conditions for suppression of FI and F2 layers as 
outlined by Chamberlain and McElroy (1966) are required* 

Since physically significant model atmospheric structures 
were adopted in these studies, it would appear that the 
E- layer interpretation of the Martian ionosphere is basically 
correct* One piece of experimental evidence, however, raises 
serious questions with this explanation* The coefficient for 
dissociative recombination of C 02 ^ ions with electrons has been 
recently measured at (3*6 t 0*5) x lO"*^ craVsec at a temperature 
of 300 ® K and appears to be of the same size at a temperature of 
^ 210® K (Weller and Biondi, 1967). This value is substantially 
smaller than the minxumim value of 10“^ cm^/sec required by the 
E-layer interpretation for suppression of other ionospheric 
features and, thus, the need for considerable modification of 


\ . 
'I' 



54 


the ionospheric models of Chamberlain and McElroy (1966) and 
KcElroy (1967) is evident* 

FI- layer interpretation , 

Donahue (1966) suggests the possibility that the ionospheric 

layer observed on 14ars is an PI region* Assuming that the 

thermal structure above 50 km obtained by Chamberlain and 

KcElroy (1966) is accxiratCi Donahue (1966) points out that, if 

the lov/er atmosphere of l^rs consists primarily of CO2 rather 

than the 44^ considered by Prabhakara and 

Hogan (1965), the smaller scale heights in the lower atmosphere 

could result in a gas density of 'V 10^^ particles/cm^ near 

125 km. In this case, he finds that the maximum photionization 

rate would occur when the sun is ^ 20® above the horizon, as 

it was during the occulta tion experiment. For a characteristic 

photoionization rate of 1,2 x 10^ cm“^3ec the observed 

electron density at the peak (9 x 10^ electrons /cm ) requires 

3 -1 

an effective recombination coefficient of 1,5 x 10 cm sec , 
This corresponds to the value anticipated for 0 ^ ions (Biondi, 
1964), and Donahue (1966) concludes that therefore 

responsible for the observed maximum* He also demonstrates 
that the suppression of an F2 peak can readily occur if the 
rate coefficient for the reaction between O'*’ and CO^ to form CO 
and 02^ is 1.2 x 10“^ cm^sec*^ (Feshenfeld, et al, , 1966), In 
that case, the O/COg abundance ratio must only remain less than 

10^ up to the altitude at which diffusive loss of O'*’ ions 
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becomes important* 
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Stewart (196?) also favors the FI interpretation. By 
applying a criterion derived by Yonezawa (1965), he has shown 
that, in the ^^rtian atmosphere, the boundary between ionos- 
pheric regions governed by attachment-type and recombination- 
type electron loss probably occurs well above the observed 
ionization maximum. Therefore, the ionosphere can be described 
by a Chapman function, with electron-ion diffusion being unim- 
portant in the vicinity of the peak, Stewart (1967) has argued 

that only slight CO2 dissociation is present and therefore, 

^ * 
that ©2 rather than 0 is the dominant ion responsible for the 

observed maximum. He finds that the temperature of the ^'artian 
exosphere varies from »v 1S0®K to 550®K depending on the point 
in the solar sunspot cycle, and that the upoer atmospheric 
temperatures for low solar heating flux are consistent with the 
I-ariner IV observations provided that the ionosphere is an FI 
layer in which molecular ions (0^*) dominate. 

The FI interpretation was not ruled out by Chamberlain and 
i'iCjilroy (1966), Since it requires less extreme conditions in 
the lower atmosphere than does the F2 model, and less extreme 
assumptions about reaction rates than does the E model, the FI 
interpretation for the ionosphere of >fe.rs seems to be the most 
acceptable. No detailed calculations, however, have been 
performed to confirm the existence of the required neutral par- 
ticle density at the ionospheric level. 
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DETER14INATI0N OF AT140SPHERIC STRUCTURE 

From the preceding discussion | the strong dependence of 
the type of ionosphere observed on Mers upon the structure of 
the lower Martian atmosphere is evident* Existing models for 
the lov/er atmosphere demonstrate the total unacceptability of 
the F2 interpretation, and favor an E-type ionosphere* The 
E-model, however, has been ruled out by the recent experimental 
determination of a critically important reaction rate* Thus, 
the FI interpretation seems to be the only plausible explana- 
tion of the observed ionosphere, but the considerable rein- 
forcement of this viewpoint, which would result from the 
actual construction of a consistent lower atmospheric model, is 
lacking. The present research is designed to determine whether 
a physically sound model for the lower Martian atmosphere can 
indeed produce the conditions required for production of an FI 
ionospheric layer at an altitude of 125 km* 

Evaluation of Model I of Prabhakara and Hogan {1% H 

In the calculations of Prabhakara and Hogan (1965), various 
assumptions were made which limit the acceptability of their 
results as realistic representations of conditions in the lower 
Martian atmosphere at the time and place of Itoriner IV occulta- 
tion* Some of these assumptions involve parameters for 
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variation, such as surface temperature, surface , ..i 

atmospheric composition, while others, such as the adoption of 
the Lorentz line shape and the neglect of non-local thermo- 
dynamic equilibrium effects, concern the basic physics employed* 
It is unlikely that any of these assumptions would drastically 
influence the results j the atmospheric structures obtained by 
Prabhakara and Hogan must be reasonable, though not necessarily 
realistic. It is possible that the cumulative effect of these 
assumptions could produce overestimates of the atmospheric 
density at the ionospheric level, and thus lead to the false 
conclusion that the Martian ionosphere is of the E-type. 

In their Model I (Figure 3), Prabhakara and Hogan (1965) 
adopted a ground temperature of 230®K, characteristic of the 
daily mean ground temperature near the I'lartian equator. 

Figure 1 (Mintz, 1961) suggests that a more representative 
value for the surface temperature at the occultation point on 
Mars (55® 3) on an early winter afternoon would be a' 180® K, 
/Although this value is based on a somewhat subjective analysis 
by Mintz (I 96 I), it is probably a better estimate of the actual 
ground temperature at the point of occultation than is the 
arbitrarily selected value used by Prabhakara and Hogan. Any 
significant reduction of surface temperature in a model atmos- 
phere v/here the structure is governed by radiative transfer 
will result in a general lowering of temperatures throughout 
the model up to great heights above the surface. 
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The surface pressure and atmospheric composition adopted 
by Prabhakara and Hogan in the construction of their Model I 
can now be replaced by more realistic estimates of these para- 
meters. The ambient surface pressure on Mars has been found 
to lie within the 4 mb to 7 mb range by the Mariner IV experi- 
menters (Kliore, ejt al., 1965) and within the 5 mb - 13 mb 
range by the spectroscopists who prefer a value of 6 mb 
(Belton and Hunten, 1966). Thus, it is likely that the atmos- 
pheric pressure at the surface of Mars is somewhat lower than 
the 10 mb value adopted by Prabhakara and Hogan. Moreover, the 
same observations show that the Martian atmosphere contains a 
much larger fraction of CO 2 than the 445^ assumed by Prabhakara 
and Hogan, with a 100^ CO^ composition being allowed as a real 
possibility (Belton and Hunten, 1966). 

Since the available evidence now indicates that the ambient 
temperature and pressure at the surface of Mars near the occulta— 
tion point were somewhat lower, and the mean molecular weight 
of the atmosphere somewhat higher than the values adopted in 
the Prabhakara and Hogan model, it follows that the actual 
scale height in the lower atmosphere was probably significantly 
smaller than the model scale height at corresponding levels. 

For a temperature of 230*K and a mean molecular weight of 
35*1 g/mole (corresponding to a 445^ CO^ “ 56% N 2 mixture), a 
model scale height at the surface of /v 14*4 km is obtained. 
Hov/ever, for a temperature of 160*K and a mean molecular weight 
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of 44 g/mole (corresponding to pure CO2), a scale height of 
A' 9 km near the surface results, in close agreement with the 
value of 6 km - 10 km measured by Mariner IV (Kliore, et al . . 
1965) for the lower regions of the Martian atmosphere. This 
reduction of scale height, coupled with a moderately smaller 
surface pressure, would result in smaller number densities 
throughout the region of the atmosphere of interest, with the 
cumulative effect of the reduced scale height producing signifi- 
cantly lower number densities in the upper layers of the model. 
Thus, it appears that updating the adopted input parameters to 
correspond more closely to the lower atmospheric conditions at 
occultation v/ould favor the FI ionospheric interpretation by 
producing a generally colder, heavier model atmosphere, with a 
decrease in atmospheric density at higher levels from the 
values previously obtained by Prabhakara and Hogan (1965), 

Another assumption, made by Prabhakara and Hogan (I965) in 
their treatment of the radiative transfer problem, probably 
results in overestimates of temperatures in all of their models. 
In their study, vibrational relaxation of CO2 was not considered, 
and local thermodynamic equilibrium (LTE) conditions were 
assumed to be present throughout the atmosphere, both in the 
lyy "6^ near infrared bands and in the 15^ far infrared band 
of CO2, A discussion of the breakdown of LTE in the vibration - 
rotation bands of CO2 is given in Appendix B, and a method 
derived by Prabhakara and Hogan (1966) for treating the radiative 
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transfer problem under non-LTE conditions is described in 
Appendix C. As shown in Appendix C, if a strong source is 
present to supply energy directly to the kinetic mode in an 
atmospheric region where the band responsible for the radiative 
structure of the atmosphere (15yk CO 2 band) undergoes relax- 
ation, large heating of the atmosphere with a consequent 
elevation of the temperature structure will result locally. 
However, as also shown in Appendix C, if no external heat source 
exists, the atmosphere will eventually revert to LTE conditions, 
regardless of the collision frequency. 

As illustrated in Figure 4, Prabhakara and Hogan (I 965 ) 
concluded that, under LTE conditions, the deposition of solar 
energy in the Martian atmosphere could raise the temperatures 
in the upper levels by as much as 20®K above the values obtained 
from basic radiative — convective considerations. Therefore, 
one might expect the temperature profiles obtained by 
Prabhakara and Hogan to considerably underestimate the temper- 
atxires in the higher Martian atmosphere where LTE breaks down 
(above /v 65 km in their Model I). However, the solar energy 
input at levels above 20 km in the models of Prabhakara 
and Kogan is due almost entirely to absorption in the 1^ - 6^ 

near infrared bands of C 02 » these bands, like the 1^ 
band, undergo relaxation at low pressures. If the solar 
energy in these bands does not pass into the kinetic mode, as 
v/as implicitly assumed in the study of Prabhakara and Hogan 
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(1965), but is scattered instead, then it cannot influence the 
temperature structure. In this case, only relatively weak 
sources and sinks of energy will remain outside of the 15yi^ 

CO 2 band, and the atmospheric temperature structure should not 
depart substantially from that obtained under LTE conditions. 
Moreover, if the near infrared energy is scattered rather than 
absorbed, the thermal structures presented by Prabhakara and 
Kogan probably constitute overestimates, rather than under- 
estimates, of the atmospheric temperatures. From these consid- 
erations, it would not be unreasonable to conclude that the 
temperature near the top of their Model I is too warm by 
^ 20*K, and that temperatures would decrease to a value of 
^ 140 ® K in this model, if the near infrared solar energy 
source were discounted. 

It should also be pointed out that the mean solar energy 
absorbed in the Prabhakara and Hogan models was calculated for 
equatorial equinox, rather than for the high-latitude, winter 
conditions encountered in the Mariner IV occultation. The 
larger mean zenith angle and shorter sunlit portion of the day 
appropriate to the occultation point on Mars imply a smaller 
solar energy input than would be present at equinox near the 
Martian equator. Modification of the solar energy input in 
the models of Prabhakara and Kogan to more closely represent 
occultation conditions would therefore have a further tendency 
to produce lower atmospheric temperatures. 


r 



62 


An additional assumption made by Prabhakara and Hogan (1965) 
which may significantly influence the temperature structures 
obtained, concerns the evaluation of their atmospheric trans- 
mission functions in the near and far infrared CO 2 bands* The 
neglect of Doppler broadening of the spectral lines results in 
underestimation of the line widths by several orders of magni- 
tude in the upper portions of their models. As a result, the 
opacity of the higher atmospheric layers has been greatly under- 
estimated, so that the temperatures in the upper portions of 
their models may have been underestimated also* 

Plass and Fivel (1953), from their study of the relative 
importance of Lorentz and Doppler broadening effects on the 
radiative transfer in the terrestrial atmosphere, concluded 
that the Doppler effect can be neglected for both strong and 
v/eak lines up to the 50 km altitude, where the ambient pressure 
is /V *9 mb, but that for lines of intermediate intensity, this 
effect may become important at lower levels* Elsasser (I960) 
suggested that the neglect of Doppler broadening for all lines 
nay be a reasonable approximation up to 'V kO km in the Earth *s 
atmosphere, where the' pressure is /v 3 mb* Young (1964) has 
found the conclusions of Plass and Fivel (1953) for the strong 
and v;eak lines to be valid; however, he claims that the impor- 
tance of lines of intermediate strength cannot bo neglected, 
and that use of the Lorentz line shape at pressures lower than 
/V 20 mb can introduce appreciable errors* 
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Following the arguments of Plass and Fivel (1953) and 
Slsasser (I960), Prabhakara and Hogan (1965) had reasoned that 
the pressure in the Martian atmosphere at which Doppler broad- 
ening becomes important is ^ 6 x 10"^ mb. In view of the 
conclusions of Young (1964) » however, it appears that this 
broadening mechanism becomes important at higher pressures 
(lower atmospheric levels). In any case, if meaningful temper- 
ature proflies, determined by radiative transfer, are to be 
obtained up to the level of the ionosphere on Mars, then the 
great increase in the opacity of the upper layers which would 
result from consideration of the Doppler broadening mechanism 
must be taken into account. 

Evaluation of the results of Prabhakara and Hogan (1965)| 
particularly their Model I, leads to the conclusion that, 
although the thermal structures they obtained are not 
unreasonable, changes in various input parameters and consider- 
' ition of the important effects of Doppler line broadening and 
vibrational relaxation of C0« would, nevertheless, produce more 
meaningful atmospheric models. Furthermore, it would seem 
that these revisions could quite conceivably create a cooler, 
more compressed atmospheric structure, with significantly 
lower densities at higher levels, and thus increase the likeli- 
hood that the ionosphere observed on Mars by I'lariner IV is of 
the FI -type. 
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Present Method 

In view of the numerous shortcomings of existing theoretical 
models, discussed in previous sections, the present study is 
undertaken to obtain a better understanding of the conditions in 
the lower ^!h.rtian atmosphere, at the time and place of I'feriner 17 
occultation. 

The present determination of the atmospheric structure on 
Kars is based upon a steady state solution of the equation of 
radiative transfer, allowing for vibrational relaxation of the 
radiating gas, and for the absorption of solar energy in the 
atmosphere. Possible relaxation of the vibration - rotation near 
infrared bands of CO 2 in which solar energy is absorbed is consi- 
dered here, and atmospheric transmission functions in the near 
and far infrared bands of CO 2 are calculated using a modified 
version of the Elsasser band model, taking both Lorentz and 
Doppler broadening of the spectral lines into account. Moreover, 
the required input data are chosen from the most recent estimates 
of the surface parameters and atmospheric composition available 
at this time. 

Radiative transfer . 

The physical theory related to the present construction of 
atmospheric models from radiative transfer considerations is 
developed in the Appendix. In Appendix A, a formal solution of 
the equation of radiative transfer is presented for boundary 
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conditions which are meaningful for studies of planetary atmos- 
pheres. In Appendix B, vibrational relaxation of the CO^ mole- 
cule, and the relevance of this phenomenon to the problem of 
atmospheric radiative transfer is discussed. In Appendix C, 
radiative transfer under LTE and non-LTE conditions is consi- 
dered, and convenient expressions are obtained for the atmos- 
pheric source function and the energy balance relationship, 
applicable v/hether or not the gas is relaxing (Prabhakara and 
Hogan, 1966). In Appendix D, the method of evaluating the 
atmospheric transmission functions required for solution of the 
radiative transfer equation is described in detail. 

As shown in Appendix A, the radiative transfer eqixation can 
be applied to studies of a plane parallel atmosphere of finite 
thickness, if appropriate boundary conditions are specified. 

The transfer equation then has a solution in terms of th ward 

and downward directed specific intensities; these can b 
evaluated at any point in the atmosphere, for any prescribed 
atmospheric structure, provided the emission of the gas can be 
related to local conditions, and the transmission of the gas 
can be calculated for an arbitrary atmospheric path. With a 
knov;ledge of the upward atnd downward intensities at every level 
and the transmission characteristics of the gas, the amount of 
atmospheric radiation absorbed per unit volume at every level 
in the atmosphere can be obtained. The atmospheric radiation 
absorbed at any level, plus the solar radiation deposited at 
that level, must, under steady state conditions, be balanced by 
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the local emission of the gas. By relating the emission of the 
gas to the temperature (taking vibrational relaxation into 
account, if necessary), a new temperature estimate at every level 
in the atmosphere can be obtained. The resulting thermal struc- 
ture will produce the emission required to balance the total ab- 
sorbed energy per unit volume, calculated from radiative transfer 
considerations for the old atmospheric structure. Repetition of 
this procedure will eventually result in a convergent solution 
for temperature as a fxmction of height in the atmosphere. 

V7ith the boundary conditions discussed in Appendix A of zero 
downward flux of infrared energy at the top of the atmosphere z^ 
and a blackbody surface at the ground Zg, the upward and downv^ard 
intensities at level z^^ are, from equations (A. 24) and (A.25)i 


iMv,G,Z3^) 


and 


B(v,Zg)^(V,e,Z^-Zg) 

/‘Zi . . i^(v,e,z^-z) 

+ / J(v,z) — ^ dz 

Jz^ cz 



^t(v,Q>z-zi) 


dz 




( 1 ) 


(2) 


The meaning of the various symbols in equations (1) and (2) is 
given in Appendix A. 

Consider a thin slab of the atmosphere of geometric thick- 
ness iz centered at the level Zq^ (Figure 9)» The energy origi- 
nating at levels below zi-(iz/2) which is absorbed by the slab 
per unit time, A^(iz), will be given by 





6a 


(^z) = yyi^(V,e,z^-(iz/2)) X (1 -r{V,0,iz)) dwd^) , (3) 

V CJ 


where (V,©,^z) is the transmission of the slab at frequency 
V and zenith ansle 6, and where the frequency integration is 
performed over the 15 ^ CO 2 band, and the angular integration is 
performed over the lower hemisphere. A similar expression for 
the downv/ard travelling radiation absorbed by the slab is 


(^z) = z-^{iz/2)) X (1 -'t(V,G,iz)) 

Ci> 


dwdv 


(4) 


where the angular integration is over the upper hemisphere. 

The total energy emitted by the slab per unit time under 
conditions of local thermodynamic equilibrium R(iz), is a 
f\mction of the mean temperature of the slab, here taken to be 
the temperature at level z^: 


R(^z) *= ^B(v ,z^) X (1 dwdv , 


(5) 


V 


where the integration with respect to solid angle is performed 


over the entire sphere. 

The mean solar energy absorbed per unit time by the slab 
Q,{£z) in the Ij^ - 6^ near infrared bands of CO 2 may be 
expressed as 


QUz) 




q>’ + X’(z) , 



V SI 


Q(V,z^) ^ ^ ^ 


v/here X’(z) and (p* are a mean relaxation time and a mean 
radiative lifetime, respectively, for the 1 ^ - 6^ bands 
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(Appendix B), Q(-v,z^) is the solar energy of frequency v incident 
per unit area per unit time at the top of the atmosphere, fre- 
quency integration is performed over the - 6^ infrared bands 
and integration over hour angle is made between noon and sunset. 
The hour angle Si and zenith angle 0 are related by the expression 

cose - sin(l) 8in(d) + cos(l) cos{d) cos St (7) 


where 1 is the latitude (55® S) and d is the solar declination 
(15.3* N). The near infrared solar flux at the top of the I^Iartian 
atmosphere Q(^),z^) was obtained from the Smithsonian Keteorolo - 
gical Tables (compiled by List, 1963) by multiplying the values 
tabulated there for Earth by the appropriate factor ( 0 . 432 ) to 
compensate for the difference in the mean distances of the Earth 
' and ^brs from the sun. 

Under steady state conditions, the total energy of solar and 
terrestrial origin absorbed by the slab must be balanced by the 
total energy emitted. From the considerations of Appendix C, we 
are able to relate the emission of the gas under non-LTE condi- 
tions to the LTE emission (equation (C. 3) ) . Also, we are able to 
represent the non-LTE effect in the 15yi^ CO 2 band as a factor 
amplifying the solar energy absorbed (equation (C.ll)). In view 
of expressions (3), (4) and (5) above, equations (C.3) and (C.ll) 
were rewritten as 



( 6 ) 

(9) 

V 
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From a simultaneous solution of (3) and (9) , the temperature 
structure of the atmosphere was obtained as follows. An arbi- 
trary initial atmospheric structxire was specified, and the atmos- 
phere was assumed to be in LTE (i.e* J(v»Z]^) “ B(v,Z]^)) through- 
out. For this structxjre, (iz), A^(iz) and Q(^z) were 
calculated from eqixations (3)» (4) and (6) and R(iz) obtained 
from (9)* Since R(Sz) is directly related to the Planck emission 
B(v,z^), which, in turn, is a unique function of the temperature 
at level z^, the height distribution of R(iz) determines a new 
distribution of temperature. The inversion of R(^z) to obtain a 
temperature was performed using a standard numerical technique, 
the Newton-Raphson method. This second temperature structure was 
then used to evaluate the source function J(v,z^) from equation 
(3), and the values of J(v,Z 3 ^) so obtained were used to solve 
equation (9) for a third atmospheric structure. This procedure 
was continued until an equilibrium thermal structure was obtained 
satisfying the convergence criterion that the temperature at any 
level varies by no more than 0.1*K in two successive iterations. 

In this manner, a temperatxire , pressure and density stratifica- 
tion of the atmosphere consistent with the source functions and 
absorbed solar energy was obtained. 

C onvection . 

In the course of these thermal calculations, whenever the 
lapse rate of temperature obtained exceeded the adiabatic lapse 
rate (g/cp, evaluated for the acceleration due to gravity on Mars 
and the chosen composition), it was replaced by the adiabatic 

c, 

r 
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lapse rate and the calculation continued* As pointed out by 
Arking (1962), this procedure for treating convection near the 
planetary surface amounts to assuming maximum efficiency for 
convective heat transfer* The restilting tenq>erature profile is 
somewhat artificial at the tropopause, since the discontinuity 
of temperature obtained there would be removed in the atmosphere 
by the penetration of convection into the stable region above. 

It is also artificial immediately adjacent to the ground where, 
in the atmosphere, the decrease in the efficiency of convective 
heat transfer with height tends to produce a superadiabatic 
gradient. Arking (1962) concludes from mixing length theory that 
the effect near the gro\ind is not serious since the temperature 
gradient approaches the adiabatic gradient a very small distance 
above the ground, so that the temperature obtained near and at 
the ground by assuming exactly adiabatic conditions should be 
correct to within a few degrees* He also indicates that if the 
region of transition from convective to radiative regimes is 
optically thin, so that the radiative flux entering the radiative 
region originates primarily below the transition region, then it 
is not necessary to knav the temperature near the tropopause 
accurately to determine the radiative structure above. If this 
transition region is optically thick, however, the situation be- 
comes more complicated. Since no means are available, in the 
absence of observations, for determining the thickness of this 
transition region in the atmosphere of ^felrs, it becomes necessary 
to ignore the effects of penetration convection and to determine 
the boundary between convective and radiative regions by the 
above procedure* 
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Acceptability of method . 

The general accuracy of this procedure for the construc- 
tion of a model atmosphere with a structure governed by 
radiative transfer in the 15yi/ CO2 band, was demonstrated 
by Prabhakara and Hogan (1966), who used this method in 
their study of the equilibrium temperature distribution in 
the Earth’s mesosphere and lower thermosphere. In Figure 10, from 
Prabhakara and Hogan, (1966), a comparison is made between 
observed summer and winter temperatxire profiles at 
Fort Churchill and White Sands, and three thermal structures 
calculated by the same basic method as outlined above. 

Prabhakara and Hogan (1966) determined profiles of temperature 
in the Earth’s atmosphere for three different values of the 
relaxation time in the 15yH CO2 band; lO**^ sec, 5 x 10"^ sec 
and 1.5 X 10”^ sec. Solar heating in the ultraviolet and 
visible by O2 and photochemically produced O3 was considered 
in their study, rather than the near infrared heating important 
in the present study of the Fartian atmosphere. As can be seen 
from Figure 10, the profiles they obtained are in good general 
agreement with the observed temperature profiles. At higher 
levels, the value of X used exerts an important influence on 
the temperature, due to the sizable amount of solar energy 
absorbed there by O2. Were this source of energy at higher 
levels v/eaker, the particular value of X chosen would have 
a smaller effect on the temperature structure. In any case. 




Temperature (*K ) 


Figure 10 « Comparison of calculated temperatures with 

observations at ft. Churchill and White Sands, 
I: X« • 10**|sec; II; » 53dO“5sec; III; 

» 1,5x10"^ sec, (After Prabhakara and 
Hogan, 1966), 
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discrepancies betv/een observed and calculated temperatures 
at higher levels can be easily accounted for by the fact 
that important physical processes such as dynamical heating 
and cooling and OH airglow emission v/ere not considered. 

The overall ability of this procedure to accurately repro- 
duce the observed temperature structure in the Earth’s 
atmosphere renders it acceptable as a method for the 
construction of model atmospheres for Mars also. 

Photochemistry . 

As is well knovm, CO2, the principal constituent of the 
Martian atmosphere, is dissociated by solar ultraviolet 
radiation into CO and 0, and, as a result of various recombi- 
nation processes, O2 and 0^ may be produced. This photo- 
chemical activity results in a decrease in the concentration 
of the constituent which is largely responsible for the thermal 
structure of the lower atmosphere. It also results in the 
deposition of solar ultraviolet and visible radiation in the 
higher atmospheric layers. Thus, it is important to assess 
the ir.fluence of CO2 photochemistry on the atmospheric structure. 

In the calculation of a thermal structure by the method 
outlined above, the composition of the atmosphere v^s assumed 
to be uniform \vith height. The possible effects of CO2 
photochemistry on this basic structure were estimated, however, 
from a separate calculation, for each thermal model, of a 
photochemical equilibrium distribution of COg a^d its 




V 
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dissociation products. Estimates of the solar heating 
produced by this photochemical activity were obtained simul- 
taneously with the distribution of the gases. Ideally, the 
structures of temperature and composition should be made 
totally consistent with each other by iteration. However, 
the distribution of gases in the Martian atmosphere is not 
determined by photochemical activity alone. Certainly, both 
eddy mixing and diffusive separation of gases in the 
gravitational field of the planet exert important effects 
on the local composition. Therefore, rigorous computation 
of consistent profiles of temperature and photochemical 
equilibrium composition is of little more than academic 
interest, and therefore such an iterative calculation was 
not performed. 

The photochemistry of CO2 in the Martian atmosphere 
has been discussed by ^^armo and Warneck (1961). The follow- 
ing eight reactions involving photodissociation and molecular 
collisional recombination were considered in their study: 


O 2 h V 

— 

0+0 

0 + 0 + M 

— 

O 2 + M 

0 + O 2 + M 


0^ + M 

0^ + hv 

— 

O 2 + 0 

0^+0 

— 

202 

COo hV 


CO + 0 

CO + 0 + M 


CO 2 ♦ M 

CO + 0 

— 

CO 2 


^1 

Kl 

a^ 

K4 


(X< 2450 A) 


(X< 11000 A) 


(X< 1300 A) 
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Here a2 and are the photodissociative production rates, 
the number of quanta absorbed per molecule per unit time by 02^ 
0^, and CO2, respectively. 

The assumption of photochemical equilibrium at every level 
in the atmosphere leads to the relations: 


“^3*^00^ * ^5^ ”cd^0 " ® 

“^l’^02 * ^2^03 2^3”03”0 “ ° 


( 10 ) 

( 11 ) 

(12) 


Furthermore, it can be seen that the total number of 0 atoms at 
any level contained in the 0, O2 and O3 forms must be equal to 
the total number of 0 atoms released by CO2 dissociation. For 
every molecule of CO2 dissociated, one 0 atom is liberated, so 
that the total number of 0 atoms released in this way must be 
equal to the number of CO molecules produced. Thus, at any time 


^co “ ’"O ^ 2no2 ^ 3no^ 

Similarly, the total number of CO2 molecules present at any level 
must be equal to the total number available nQQ^(o), less the 
number of CO molecules produced by dissociation: 


^C02 “ *^002^°^ " "CO 


(14) 


Equations (10), (11), (12), (13) and (14) represent five 

relationships between five unknown quantities, nQ, nco 

and n ^ so that a simultaneous solution for these quantities 
CO2* 


c. 

V 
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IS possible. By a succession of substitutions and rearrangements 
an expression for the number density of atomic oxygen n^ in terns 
of known quantities can be derived: 

^ 0 ^ ^ AnQ^ -I- Bno^ + Ouq + D « 0 (15) 

where 

^ ^ 2KiK 3« 3KiK2nM) - (K4nM-«- K5)(aiK3-»- aag^n^^) 

E 

g » ~^3 ^ ^2^3^00? ^ aiK 3 ) - ( + K 5 )aia 2 

E 

Q „ ^3^^A^C02(o) 

E 

E 

and 

E - i1k( \n^ ♦ Kj ) ( KjKj - 2Kj^Kj - 3Kj^K2n„) . 

This quartio algebraic equation may be solved analytically or 
numerically for n^. 

V/ith a knowledge of Uq, the O 3 number density can be 
evaliaated from 

KiKjnc^n^ 

( 16 ) 

^ ai(a 2 + K3no) - K2K3nQ2nM 

and, with a knowledge of Uq and the O 2 number density can 
be calculated from 


c. 

y 



(17) 


n 


O 2 


(a^ + K^nQ)nQ^ 


Finally, the number densities of CO and CO 2 can be obtained from 
equations (13) and (14). 

The photodissociation rates (aj, ’s) required for the evalu- 
ation of the above quantities may be expressed at any level as 


(zi) = J (^i(V)F(v,z.p)exp^sec6 ^ Oj(V)y^ 


n^ ( z)d^^dV 


(IS) 


vjhere CTj^(v) is the dissociation cross section of the i consti- 
tuent at frequency V, and F(V,Ziji) is the incident photon flu>c per 
unit frequency at the top of the ^tLrtian atmospliere. Frequency 


integration is carried out over the range in which the dissoci- 
ation cross section of the i^^ constituent has non-zero values. 

The mean solar energy absorbed in the layer by the 
dissociation of CO 2 , C '2 and 0^ may be expressed as 



Q’(iz) = ^ //Q»(V,Zt) 


V SI 


3 

exp (-sec © CT-j (v) 


<- 

/- 




n^ (z)d: 


-exp(-sec0 "511 O’-; (v ) 


z-j_-(4z/2 ) 


n -; ( z ) d z 
J / 


(19) 

dsidv 


vihere Q»(v,z^) is the solar energy of frequency V incident at 
the top of the atmosphere, frequency integration is performed 
over the relevant ultraviolet and visible spectral regions, and 
integration over hour angle 51 is made between noon and sunset, 
as in equation (6). The relationship between the energy flux 
and photon flxrx at the top of the atmosphere is given by 
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In the present photochemical calculations, the spectral 
distribution of solar energy was obtained from several sources* 
At wavelengths shorter than 1800 A an energy distribution from 
Hinteregger, ^ al. (1965) was adopted, while between 1800 A and 

O 

2600 A, values of the soQ^ar flux were taken from Detwiler, et al . 
(1961), For v/avelengths between 2600 I and 7000 X, the solar 
energies were obtained from the Handbook of Geophysics (compiled 
by Cast, 1957). The values of the solar flux given by these 
sources apply at the distance of the Earth from the sun. To 
compensate for the greater I'fers - sun distance, they v;ere multi- 
plied by the appropriate reduction factor. 

The dissociation cross sections 0]|^ for O2, Cr^ for 0^, and 
0*3 for CO2, are also required for these calculations. For 
wavelengths less than I7OO A, and O3 were taken from Kicolet 
(i960). Values of between 1700 A and 2200 A were taken from 
V/atanabe (1956) and for wavelengths greater than 2200 A from 
London, ^ al, (1962), Estimates of 0*2 between 2300 A and 3OOO A 
and in the Chappuis bands in the visible were obtained from 
Vigroux (1952) while for wavelengths shorter than 2300 A, the 
values of Inn and Tanaka (1953) were used. 

The reaction rates K|j_ and K2 were taken from Ka\ifman and 
Kelso (1961) and was taken from Campbell and Nudelman (I960). 

A value of suggested by Bates and Witherspoon (1952) was used, 
and was assumed to be zero. The values adopted for the 
reaction rates are 
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= 3.0 X 10"^^ cm^mol“^sec'’^ 

K2 “ 1.2 X 10“^^ cm^mol“^sec"^ 

K 3 = 9.0 X 10‘^^ exp(-2133/T) cm^mor^sec“^ 

K|^ = 5.0 X 10’^^ cm^mol"^sec“^ 

Using these values of solar flxix, dissociation cross 
sections and reaction rates, height distributions of CO2, CO, 0, 

O2 and O3 were calculated for each of the thermal models of the 
lower atmosphere as follov/s. For the equilibrium temperature 
structure given by a particular model, the photodissociation 
rates (a^j^’s) in the uppermost layer of the atmosphere were evalu- 
ated from equation (16). Next, equation (15) v/as solved numeri- 
cally for Hq, using the Nev/ton-Raphson method of converp^ence , and 
the number densities nQ^, nQQ and ^002 determined from 

equations (16), (17), (13) and (14). These number densities were 
then used to determine new values of the dissociation rates for 
the top layer, and the cycle was continued until a convergent 
solution for the concentration of each gas in this layer v/as 
obtained. The solar energy absorbed in the uppermost layer was 
then calculated from equation (19), for this photochemical equi- 
librium distribution of gases. The calculation then proceeded to 
the layer immediately below, and the entire sequence was repeated. 

Height distributions of composition and short-wave solar 
heating, consistent with the temperature structures calculated 
earlier, were obtained in this manner for each model. The 
validity of the temperature structures was examined in view of 
these photochemical considerations. The estimated effects of 
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photochemistry, eddy mixing and gravitational separation on the 
atmospheric models is discussed in the presentation of the 
results. 


Numerical procedure . 

In the evaluation of the mean rate of solar energy 
absorption in a given atmospheric layer Q(v>z) and i^’(V,z) from 
equations (6) and (19), integration over hour angle from noon 
to sunset v;as approximated by a summation over nine equal hour 
angle intervals, each v/ith a characteristic zenith angle. The 
sunset hour angle was calculated from equation (7). At 
sunset the zenith angle 0 is 'ir/2, so that, from (7) 


» arc cos ^-tan(l) tan (d)^ 


( 21 ) 


For the Mariner coordinates SI is /v 67* so that each hour 

s 

angle interval is ^ 7,5*. The factor l/SV v/hich precedes the 
integrals in these equations is required to obtain the mean 
value of the solar energy input over the total day ( 5t “ 2TI') . 

In the evaluation of absorption and emission in the 15^ 
CO2 bands from equations (3), (4) and (5), cylindrical syimaetry 
about the vertical was assumed. Integration over solid angle 
may then be replaced by integration over zenith angle, if the 
factor dfl on the right hand side of these equations is replaced 
by the factor 2TTsin©d©. ©-integration is performed from 1T/2 
to iy in equation (3), from 0 to tf/2 in equation (4) and from 0 
to Tr/2 in equation (5)» where an additional factor of 2 must be 
included to take into account the fact that the layer radiates 
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over a sphere. In these integrations, nine angular steps oi 10“ 
each were used. 

To test the sensitivity of the results to the height 
resolution of the models, calculations v/ere performed using 
various layer thicknesses, ranging from 1 km to 10 km, i.lco, 
in some cases, the layer thickness was varied with altitude, 
with very thin layers (0,1 km) adjacent to the ground, i.o 
significant differences in the calculated temperature struc- 
tures were observed as the height grid was varied within the 
0,1 km to 10 km range. However, it was found that, since ;.iOst 
of the numerical computation is involved with the evaluation of 
the required arrays of far infrared transmission functions, the 
time required by the coi^puter to establish a convergent temper- 
ature structure increases approximately as the square of the 
number of layers. Use of layers of less than 1 Ion thickness 
v/as found to be uneconomical, while use of very thicK l;iyors 
resulted in the loss of interesting aetail (such as the convec- 
tive layer), A compromise was therefore made, v.'ith layers of 
2 km, 3 krn and 4 km thickness being used. Most of the models 
presented in the results (Models I thru VII) were calculated 
with a height scheme consisting of fifteen 2 km layers between 
the surface and 30 km altitude, ten 3 km layers between 30 km 
and 60 km altitudes, and ten 4 l«n layers between 60 km and 100 km 
altitudes. Some, hov/ever, (Models VIII and IX) were obtained 
for sixty uniform layers of 2 km thickness between the surface 
and the 120 km level. 
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An investigation was also undertaken to determine the 
effect on the calculated temperature distributions of refine- 
ment of the spectral intervals in the region. The 

strongest portion of this region (700 cm"^ to 600 cm”^) v/as 
divided first into two 50 cm ^ wavenumber intervals, then into 
five 20 cm**^ wavenumber intervals. For each case an equilibrium 
structure was obtained for the lower 100 km of the Earth’s 
atmosphere. The two temperature profiles v/hich resulted dif- 
fered by, at most, 3*K within the height range of the model, 
and both were in good agreement with normally observed tempera- 
ture distributions (Figure 10), It was therefore deemed 
acceptable to use the broader spectral intervals in the present 
thermal calculations for the Martian atmosphere, in the interest 
of economy of computer time. In Table D,I of Appendix D, two 
groups of band parameters derived from the transmission tables 
of Stull, e^ al, ( 1963 ) are presented for the 15yq 00 2 bands. 

The second group, corresponding to broader spectral intervals, 
was adopted for the calculation of the necessary transmission 
functions in this study. 
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Lower atmosphere . 

In all, nine basic models of the lower Martian atmosphere 
were constructed. The values of surface temperature and pres- 
sure and atmospheric composition adopted in each of these 
models are listed in Table III, Surface temperatures within 
the range 150*K to 220*K and surface pressures within the range 
4 mb to 6 mb were considered. In five of the models, the atmos- 
phere v/as assumed to be composed entirely of C02> while in the 
remaining four, an ^Qfo CO^ - 20;^ N2 mixture was adopted. 

The height distributions of temperature obtained for these 
cases are shovm in Figures 11, 12 and 13 • In these figures, 
the computed temperature structures are compared with those 
given previously by Johnson (1965) and Prabhakara and Hogan 
(1965), The temperature structures obtained for Models VIII 
and IX v/ere found to be nearly identical (within 1“K at all 
levels), and in Figure 13, only one curve, that of Model VIII 
is shovm to represent both models. 

It is immediately apparent from Figures 11, 12 and 13, 
that while all of the model atmospheres obtained here are 
significantly warmer than the model of Johnson (1965), they 
are significantly colder than Model I of Prabhakara and Hogan 
(1965), Only Model IV of the present study attains a tempera- 
ture slightly in excess of the Prabhakara and Hogan model 




TABLE III 


BASIC CHARACTERISTICS OF MODEL I-MRTI/J'; A'M-.OSPHE...^, 


Model 

Siirface 

Temperature 

Surface 

Pressure 

Surface 

Number Density 

Conpo 3X 



("K) 

(mb) 

(10^*^ mol/crP) 

(VJGO 2 ) 

( 2 

I 

160 

4 

l.Sl 

100 

0 

II 

ISO 

5 

2,01 

100 

0 

III 

200 

6 

2,17 

C) 

0 

c 

IV 

220 

7 

2,31 

100 


V 

150 

4 

1,93 

so 

20 

VI 

175 

5 

2,07 

so 

20 

VII 

200 

6,3 

2,26 

so 

20 

VIII 

ISO 

S 

3.22 

100 

0 

IX 

ISO 

S 

3»22 

so 

20 
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Ficure 11. Fiartian lower atmosphere: Models I thru IV 
— ^ (lOOJ'o C0«). J: Johnson (1965); P'ScH: 

Prabhakara and Hogan (1965) • 
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Martian lower atmosphere: Models VIII 
{loojo COo), IX (80;^ COo - 20;^ N?), and 
IX A. J: Johnson (1965); Prabhakara 

and Hogan (1965)* 
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(by /V 5*K), and this occurs only within a limited altitude 
span* On the other hand^ the present models are colder than 
the Johnson (1965) model only near the surface, with much 
warmer temperatures (by at least 40*K) in the upper layers. 

Models I, II, III, IV and VIII were constructed for an 
atmospheric composition of pure CO^, while Models V, VI, VII 
and IX were based upon an CO 2 - 20^5 mixture. By 

examining these nine basic models, as depicted in Figures 11, 

12 and 13, the effects of changes in the input parameters 
(surface temperature, surface pressure and composition) on the 
atmospheric structure can be estimated. 

From both Figures 11 and 12, it can be seen that an 
increase in surface temperature results in a general elevation 
of atmospheric temperatures throughout the model. Also, this 
increase in atmospheric temperatures with the temperature of 
the ground is less pronounced in the higher layers of the atmos- 
phere. In Figure 11, a 60®K increase in surface temperature 
flrom 160*'K (Model I) to 220*’K (Model IV) produces a rise of less 
than 30®K in the temperature at the 100 km level (131,6*K to 
15 s. 0*K). In Figure 12, a rise of 50*K in the ground tempera- 
ture from 150*K (Model V) to 200*K (Model VII) is accompanied 
by an increase of less than 25*K near 100 km (127, 9*K to 
149 . 2 *K). Thus, the range of temperature in the upper layers 
of both of these groups of models is less than one-half of the 
range in ground temperature. 



90 


It can also be seen from Figures 11 and 12, that the 
tropopause height increases in a linear fashion with the 
temperature of the ground. For the coldest value of surface 
temperature considered, 150° K (Model V, Figure 12) the tropo- 
pause level is found at 0,5 km, v;hile for the warmest 
surface temperature, 220°K (Model IV, Figure 11), the tropo- 
pause is located at ^ 4.5 km. Also, if a reasonable 
extrapolation is allov/ed, the calculations suggest that for 
a siarface temperature below a' 140° K, no convective region 
would appear. The results seem to indicate, therefore, that 
the height of the tropopause on I^ars (h^) varies approximately 
according to the relationship 

h » 4.5 (T - 140) / SO 

5 

where is the surface temperature in ° K and h^ is expressed 
in kilometers. The tropopause heights sha'/n for the nine 
models follov/ this relationship fairly well, falling in every 
case vri-thin 0,25 km of the predicted level. It should be 
pointed out, hov/ever, that, due to the use of rather thick 
(2 km) height intervals in the construction of these models, 
the exact location of the tropopause is somev/hat uncertain, and, 
therefore, the height of the tropopause shown in the figures 
represents to some degree an artistic interpolation on the part 
of the author. Nevertheless, in all the cases considered here, 
the tropopause level appears far below the 14 km height 

suggested by Johnson (1965). As can be inferred from the 
figures, radiative transfer is able to rapidly stabilize the 
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lower layers of the CO^-rich I-^rtian atmosphere, and as a 
result, only a thin convective region can develop. From 
radiative transfer, v/hich was not considered in his study, 
the thick convective layer proposed by Johnson appears 
unjustifiable. 

Comparison of Models II and VIII dem.onstrates the effect 
of variation of surface pressure from 5 mb to C mb when 
surface temperature (1S0°K) and composition (pure CO 2 ) are held 
fixed. An increase in surface pressure by 3 mb (60j) causes a 
slight reduction in temperature throughout the range of the 
models, with a maximum decrease of near the 100 km level. 

The similarity betv;een Models VIII and Jli shov/s that the 
variation of CO^ concentration from to IOQ 7 for a fixed 
surface temperature (1S0°K) and surface prossxure (S mb) has a 
very small influence on the thermal structure. 

The real significance of these temperature profiles, 
however, lies in the number densities they produce at hierher 
levels. The temperature, pressure and number density/ (on the 
basis of no photodissociation) obtained at the 100 km level 
in each model is given in Table IV, v;ith corresponding values 
from Johnson (1965) and Prabhakara and Hogan (1965). As 
indicated in Table IV, the smallest number density obtained in 
any of the models at 100 km is larger by a factor of 3 4 

than that obtained by Johnson. Also, the largest number density 
obtained at 100 km is more than an order of magnitude smaller 
than the relevant model of Prabhakara and Hogan indicates. 
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TABLE IV 


TEI^ERATURE, PRESSURE AND NTOIBSR DENSITY AT 100 KK 


Model 

T 

100 

(“K) 

^100 

(mb) 

n 

100 3 
(mol/cm^) 

I 

131.6 

2.5x10“^ 

1.4x10^^ 

II 

141.3 

9.0x10"^ 

4.6x10^^ 

III 

147.7 

2.0x10“^ 

9.9x10^^ 

IV 

156.0 

5.5x10“5 

12 

2.5x10^ 

V 

127.9 

4.4xl0~^ 

2.5x10^^ 

VI 

137.0 

1.0x10“^ 

5.5x10^^ 

VII 

149.2 

2.9x10"^ 

12 

1.4x10 

VIII 

133.4 

l.lxlO"^ 

5.7x10^^ 

EC 

134.0 

2.0x10“5 

l.EclO^^ 


Prabhakara 
and Honan 
(1965) 

156.5 

1.2x10“^ 

5,3x10^^ 

J ohns on 

65 

^ 5x10“'^ 

/V 2pcl0^° 


(1965) 


>:< 


On the basis of a homogeneous atmosphere. 


In the present models for a homogeneous atmosphere, a range in 

number density of 1.4 x 10^^ mol/cm^ to 2.5 x 10 ^^ 

was obtained at the 100 km level, as compared to a value of 
XO 3 

^ 4 X 10 mol/ cm given by Johnson and a value of 5.3 x 10^^ 

mol/cm^ given by Prabhakara and Hogan. In T'odels VIII and IX, 
calculations Mere carried to the 120 km levels, with number 
densities of 2.9 x lO^O mol/cm^ and 6.5 x 10^° mol/cm^ being 
obtained at that level, respectively# 

For a fixed composition the number density in the upper 
layers increases as the surface temperature is increased 
(I-.odels I thru IV , V thru VII). Increasing surface pressure, 
with ground temperature and composition held fixed (Kodels II 
and VIII), has little effect on the number densities at higher 
levels, since the larger density at the ground is offset by 
a faster decrease in density with height in the cooler atmos- 
phere produced. An increase in the CO^ concentration, with 
surface pressure and temperature held fixed (liodels VIII and 
IX) produces lower number densities in the upper layers, 
although the temperature distribution remains essentially the 
same, due to the decrease in scale height with the larger 
molecular weight. 

An interesting feature appearing in all of the models 
examined is the blunt temperatiire maximum in the 65 km - 75 km 
region of the atmosphere. Such a maximum might be expected 
to result from the absorption of solar energy at these levels. 
This interpretation is not correct, however, as can be seen 
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from r'igure 14. Model II A in this figure shows the thermal 
structure obtained for the same basic conditions adopted in 
Model II, in the absence of solar energy absorption in the 
atmosphere. Since a temperature maximum is obtained near 65 km 
in this model also, it must be a feature produced in some v/ay 
by radiative transfer in the 15 /f GO^ bands. The possibility 
of obtaining positive gradients of temperature under radiative 
equilibrium conditions in a non-grey atmosphere in the absence 
of solar heating has been noted by Goody {I 964 ). 

Analysis of the energetics involved in the determination 
of the atmospheric structure at these levels reveals the cause 
of this temperature maximum. The total 15yq energy absorbed 
by each atmospheric layer per unit time v;as divided into tv:o 
classifications: indirect radiation, originating in other 
atmospheric layers, and direct radiation, emanating from the 
planetary surface, iach of these classifications was further 
subdivided for examination according to spectral intervals. 
Thus, the energies of ’’atmospheric” and ’’surface” origin 
absorbed by each atmospheric layer in each of the chosen spec- 
tral intervals were obtained. It was found that an important 
contribution to the total 15 ^ energy absorbed per unit time 
in the atmospheric layers between 60 km and 70 km is the 
absorption of radiation directly from the surface in the v/eak 
intervals of the 15 ^ CO 2 region. 

Energy balance in the 15yt{ region requires that the cotal 
15 u energy absorbed per layer be balanced by the total 15 
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emission. Ordinarily, the largest contributions to both 
absorption and emission of energy in an atmospheric layer would 
be from the strongest spectral intervals near the center of the 
15 ^ region. However, an increase, by some means, of the rela- 
tive contribution of the v;eak intervals at the edges of the 15yf 
region to the total absorption of 15y< energy must be accompanied, 
under balance conditions, by a general increase in the 15^ emis- 
sion, most of which v/ill be accomplished by more intense radia- 
tive cooling in the strong central portion of this region. 

Thus, a rise in temperature must occur. 

At levels in the present models, where the relative contri- 
bution of surface energy in the weak spectral intervals increases 
with height (between 20 km and 60 km in l-.odel II A) , the total 

15 energy absorbed per layer decreases less rapidly with 
/ 

height than does the opacity in the strong intervals near the 
center of the 15yt< region, ^it higher levels, v;here the relative 
importance of the surface energy in the weak intervaxs decreases 
again (above /v 65 km in I.odel II A), the total energy 

absorbed per layer decreases more rapidly with heigut tran does 
the layer opacity near the band center, as a result, a temper- 
ature maximum is produced (betv/een 60 km and 65 km in I.odel II A) , 
The absorption of a relatively large amount of surface 
radiation in the weaker portions of the 15y^ region in the 
higher layers of these models is a result of the Doppler broad- 
ening machanism. In the Martian atmosphere where very low 
pressures are encountered, the fractional absorption of 
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individual layers lies along the Doppler plateau (represented 
by the lowest line in the network of Figure D.VIII, Appendix D) 
throughout a considerable height range, v/ithin this range, 
only snail decreases in layer opacity with height occur, since 
fractional absorption varies slov;ly with pressure and path when 
^ne center of the Doppler broadened lines becones saturated, 

3y examining tne fractional absorption per layer as a function 
Ox altitude, it can be seen that the Doppler plateau is reached 
at quite different levels by different spectral intervals, due 
to the large variation of mean line intensity from one specural 
intervai to another (Table D.I), v;ith weak intervals attaining 
Dopp-i.er sai^uration at much lower levels than strong intervals, 
nS u. consequence oi tnis Joppxer ;f.c'ect, the overall decrease 
iciyer opacity witn n-;ight is more rapid in trie lowei' Titmos- 
phcric j.ayers lor strong intervals tiian for weak intervale. 
Since the opacity in the weak intervals increases relative to 
that of the strong intervals with height, and since tlie freely 
transmitted radiative flux from the surface in the v/eak inter- 
vals is larger than the strong-interval flux of atmospheric 
origin, the importance of weak-interval surface energy absorp- 
tion increases with height relative to that of the st 2 '‘ong- 
interval atmospheric energy absorption. In the vicinity of 
60 km - 70 km in these models, absorbed energies of both types 
are of about equal weight. 

Thus, the temperature maxima appearing in these models are 
a consequence of the Doppler broadening mechanism and the 
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non-;^rey representation of the 15 ^ CO^ bands used in the present 
study. Such temperature maxima did not appear in the models of 
Prabhakara and Hogan (1965) since these authors did not consider 
the Doppler broadening effects. It can be seen from Figures 11 
and 12 that the maxima appears at higher levels and becomes 
more pronounced as surface temperature is increased. The upward 
shift is a consequence of the rise in density at these levels 
which accompanies the increase in atmospheric temperatures, 
while tlie greater amplitude is due to the larger weak-interval 
radiative fluxes from the warmer surface. 

This discussion is not intended to demonstrate that this 
temperature maximum is a real feature of the Martian atmosphere, 
but rather that a tendency toward such a maxinum exists when 
non-grey radiative transfer by Doppler brcaaened lines is 
treated. For the particular band moael and input parameters 
adopted in the present study, this temperature maximur.i is pro- 
duced. Refinement of the spectral intervals, or a different 
choice of bana parameters might remove this feature, but tne 
latent tendency toward ins appearance would still be present, 

V/ere t\\Q relative contribution of surface emission in the weak 
intervals of the 15 H region to the total absorbed 15^ energy 
decreased by some means, the general radiative structure of 
the atmosphe.-e would nevertheless be preserved, witn a lov/ering 
of atmospheric temperatures in the upper portions oi all these 
models, by ^ 5*K in the colder cases to ^ 10*K in the warmer 
ones, A small reduction in the number densities at higher levels 


would follow. 
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In all of the nine basic models, a value of 10”^ sec at a 
pressure of one atmosphere was adopted forX, the vibrational 
relaxation time in the 15y^ CO^ band. The value of <p, the 
radiative lifetime in this band is a. 4.12 x 10"^ sec (Appendix 
B). Moreover, in these basic models, the 1 yCf - near infra- 
red CO^ bands where solar energy is absorbed were assumed to 
relax in the same manner as the 15yt< region. A mean relaxation 
time X and mean radiative lifetime (p were defined for the near 
infrared bands (equation 6) and, in the nine basic models, the 
ratio X/<p v/as assumed to be the same as the ratio X/cp in the 
I5y^ fundamental. As can be seen from equations (6) and (9), 
when near and far infrared bands relax in the same way (i.e. 
y/<p = X/<p) , the solar energy source term in the balance 
relationship (9) takes the same form as if no relaxation \>rere 
present; all of the involved solar energy may be thought of as 
being absorbed and the atmosphere may be thought of as radiating 
with the Planck intensity. Since, as noted in Appendix B, it is 
likely that the \j{ - 6^ bands relax at higher pressures than 
the 15yt^ band (i.e. X/cp < X/<p ), the assumption that X/<p = X/(p» 
probably results in an overestimate of the importance of solar 
energy in determining the temperature structure. 

The influence of the absorption of solar energy on the 
models can be seen from a comparison of Models II and IIA 
(Figure 14). As mentioned above, an alternate model v/as con- 
structed for the same conditions as Model II, except that no 

solar energy input into the atmosphere was considered. The 
resulting Model II A is found to be colder throughout, but by a 
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very small amount, with a maximum decrease in temperature of 
/V 4®K in the upper layers. This slight decrease in atmospheric 
temperatures when solar energy is removed altogether implies 
that, even if near infrared scattering becomes important at 
levels far below that altitude at which the 15yt< band relaxes, 
the atmospheric temperature profile for a given model could not 
fall below this ”no solar energy" limit. 

Also shown in Figure 14 is the temperature profile which 
would result for the same basic conditions as were adopted in 
Model II if, hov/ever, solar energy is not scattered, and only 
the 15 ^ band relaxes, with a lifetime of 10“^ seconds. As can 
be seen from Model II B, under these circumstances, above 55 km, 
where relaxation in the 15 y‘< band occurs, temperatures rise to 
a value in excess of 200° K near the 100 km level. V.'ere a value 
larger than 10"^ sec used for the 15 relaxation time, non-LTE 
conditions would be reached at a lower altitude and even larger 
temperatures in the upper layers would result. Since relaxation 
certainly occurs in the near infrared bands. Model II B is of 
limited usefulness to the present study. However, it might be 
used for purposes of comparison to measure the effects of an 
external heat source, supplying energy directly to the kinetic 
mode, on the atmospheric structure under non-LTE conditions. 

A case was also considered in which the near infrared bands 
were assumed to relax at higher levels than the 15^ band (i.e. 

) • For the same conditions as were adopted in 

-5 

Model IX, Model IX A was constructed, using a value of 10 sec 
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f or X and a value of 10 ^ sec for\* As indlca'ted by Figure 13 » 
the thermal structure in this model becomes warmer than the 
Kodel DC structure in the upper layers , with temperatures rising 
from the Model IX values above ^ 40 km toward a maximum differ- 
ential of 'v io<»K (144.3‘’K vs. 134.0«K) at the top of the model. 
The ratio (1+V9 )/(!■'■ ) approaches a constant value of 10 

at higher levels in Model Aj as opposed to a value of unity 
in Model XX. Thus, more weight (lOx) is given to the solar 
energy absorption in Model IX A and the balance temperature is 
increased accordingly. The solar energy source term, in equation 
(9), Q(iz) (1+ 'X(z^)/(|)) , although xoniformly larger in the upper 
layers of Model IX A than at corresponding levels in Model IX, 
nevertheless decreases with height in the same manner as does 
the term representing 15^ absorption (A^Sz) + A (6z)). 

(Both terms are of approximately the same size in Model IX A, 
but the former is lOx smaller than the latter in Model DC.) 
Therefore, the temperature does not increase markedly with height 
in Model IX A but approaches a fairly constant value, slightly 
in excess of the Model IX temperatures. Because of the warmer 
temperatures, a slight increase in density at the 100 km level 
from 1.1 X 10^^ mol/cm^ in Model DC to 1.5 x 10^ mol/cm^ in 
Model DC A results. 

The relative values of the relaxation times adopted for 
Model DC A place "undue emphasis dn the importance of solar 
energy absorption. Nevertheless, no great increase in tengsera* 
ture due to this effect is observed. It should be emphasized 
that, since the relaxation time is probably longer in the near 
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than in the far infrared bands, the profiles obtained for the 
nine basic models are probably accurate representations of the 
thermal structure for the sxirface pressxire, surface temperature 
and composition assumed in each case, regardless of the exact 
value of the relaxation time in the 15 band. 

For each of the nine basic thermal models obtained, a 
photochemical equilibrium distribution of the atmospheric gases 
vras computed, using the values of the solar flxnc, photodissoci- 
ation cross-sections and recombination rates referred to 
previously* A consistent distribution of the solar energy 
absorbed in the photodissociation processes involved was calcu- 
lated simultaneously* In all nine cases, almost identical 
distributions of gases were obtained, with, however, slight 
upv/ard or downward displacements in the level of maximum photo- 
chemical activity, due to vertical shifts in density levels 
caused by the variation in atmospheric temperat\ires and assumed 
composition from model to model* Due to the extreme similarity 
in all of the resulting chemical structures, only one is shown 
(Figure 15 )i that corresponding to thermal Model IX. 

As illustrated in Figure 15, ^^ 2 * principal constituent 
in the lov^er atmosphere is highly dissociated above ^ 75 km 
where the region of maximum photochemical activity appears* 

The abundances of CO and 0 at higher levels are equal, and both 
decrease with altitude with the scale height of COg* As a 
result of a three-body 0 recombination process, a layer of 0^ 
is formed, with a maximum concentration of ^ 4 x 10^^ mol/cm^ 
near 75 km* The importance of the formation of O 2 from the 



Figure 15 • Photochemical equilibrium distri- 
bution of gases corresponding to 
thermal model IX. 
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dissociation of CO^ in a planetary atmosphere v/as first stressed 
by Chamberlain (1962). Since 0^ absorbs strongly in the same 
region of the ultraviolet where CO 2 absorption occurs, the 
presence of an O 2 layer such as that in Figure 15 serves to 
shield the CO^ at lower levels from dissociating solar energy. 

This effect is seen in the large decrease in the concentration 
of C 0 „ between 75 km and SO km, the region in which the 0 opac- 

O 

ity belo\'; 1S80 A becomes quite large. The total amount of O 2 
contained in this layer is /v 3*6 x 10 ^^ mol/cm^ or ^ 14 cm-atm, 
well below the 70 cm-atm upper limit suggested by the observa- 
tions of Kaplan, Munch and Spinrad (1964). 

Another three-body recombination reaction of 0^ and 0 
results in the production of a subdued 0 ^ layer which also peaks 
in the vicinity of 75 km in this model, with a maximum concen- 
tration of /V 4 X 10^ mol/cm^. The 0 /O ratio near the maxima 

^ ^ -6 

in the concentrations of these gases is thus 'v 10 , approxi- 

mately the same as the value of this ratio near the 0^ maximum 
in the Earth’s atmosphere. The total amount of 0^ present is 
only 'V 3.4 X 10^^ mol/cm^ or 1.25 x 10 ^ cm-atm. The influence 
of this 0^ layer on the concentrations of the other gases present 
and on the energy balance of the atmosphere is therefore 
negligible. 

The high degree of CO^ dissociation in the upper layers of 
this model is a direct result of the choice of K^, the rate for 
the three-body recombination of CO with 0. The ratio of CO 2 to 

O 2 in the upper atmospheric layers is proportional to 
(K,a,/ K,a,). Although a , a and K are known within narrow 

4113 13 1 
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limits, the reaction rate is not well known. Barth (1964) 

has suggested a value for of /v 1.4 x 10“^^ cm^sec, larger 

by more than two orders of magnitude than the value used in the 

present study. If Barth’s value is correct, the CO 2 profile in 

Figure 15 would be shifted toward higher number densities by 

more than two orders of magnitude, so that the concentrations of 

CO 2 and O 2 in the upper layers would be of about the same order. 

^^2 would nevertheless be strongly dissociated, since even the 

much larger value of K given by Barth (1964) suggests that the 

4 

CO molecule is stable against oxidation. 

Evaluation of the solar ultraviolet energy absorbed per 
layer in establishing the distribution of gases shown in Figure 
15 indicates that above the 75 km level, this energy source, if 
100 ^ efficient, would be much larger than the 1 ^ - 6yt{ infrared 
energy source corresponding to the same distribution, Below^ 

75 km the near infrared energy source would strongly predominate, 
hov/ever. Since the concentration of CO^ above ^ 75 km is 
sharply reduced by photodissociation, which results simultane- 
ously in the deposition of a sizable amoimt of solar energy in 
the upper atmosphere, a tempera txire increase in the layers 
immediately above the region of maximum photochemical activity 
near 75 km would be expected from considerations of photochemical 
equilibrium and radiative transfer alone. However, if the photo- 
chemical equilibrium distribution of gases sho^m in Figure 15 
were correct, with CO^ highly dissociated, then other mechanisms, 
the ground state transition 0(^p ) 0(^P ) in atomic oxygen, 

JU M 
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radiation froni the first vibrational and rotational levels of CO 
and the heat conduction process would assume the role of CO^ 
radiative transfer in establishing a balance v;ith the absorbed 
solar energy and determining the atmospheric temperature struc- 
ture. These three processes, all of a thermostatic nature, would 
prevent the buildup of any large temperature gradients above 
/V 75 km, so that even if photochemical equilibrium were a 
reasonable assumption, the temperature profiles given for Models 
I thru IX (Figures 11, 12 and 13) v/ould remain basically correct, 
with the possibility of a slight vjarming in the upper layers. 

The assumption that the atmospheric gases on I^jars are 
present in a photochemical equilibrium distribution is highly 
questionable, hovrever. The photodissociation lifetime of a CO^ 
molecule in the upper atmosphere of Mars is of the order of one 
month, while the corresponding lifetime for 0^ is of the order 
of several days. One is forced to conclude, therefore, that 
turbulent mixing and/or molecular diffusion will greatly influ- 
ence the distribution of gases. Stewart (1967), for a photo- 
chemical model very similar to that shovm in Figure 15, has 
compared the photochemical lifetimes of the various constituents 
with their diffusion times. The diffusion time is the time 
required for the molecules of a particular species to diffuse 
through a distance equal to the scale height of the constituent. 
Stewart finds that the photochemical lifetimes for CO^ and CO 
are greater than the diffusion time nearly everywhere, so that 
the photochemical equilibrium number densities of these 
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constituents are probably greatly in error. Also, at levels 
above and below the region of maximum photochemical activity, 
the profile of 0 is governed by diffusion, while the 0 ^ concen- 
tration is generally determined by photochemical activity 
except at levels below the maximum. Any removal of O 2 by dovm- 
ward diffusion, how^ever, will result in a lowering of the level 
at which CO^ becomes strongly dissociated, Stev;art points out 
that, as a result of simple consideration of the interaction of 
photochemical and diffusive processes, and acceptance of the 
value of the reaction rate K given by Barth (1964), it is 
necessary to conclude that CO 2 is dissociated to a point at 
which the total CO^ amount remaining is less than the amount 
observed, A possible explanation of this discrepancy between 
theory and observation is based upon the experimental evidence 
of Harteck, ^ al, (1966), who suggest that CO may be rapidly 
oxidized in the presence of trace quantities of hydrogen com- 
pounds, If the water vapor observed on 1‘Cars is responsible for 
such oxidation of CO, then turbulent mixing must be present up 
to the levels where CO is produced, so that CO can gain access 
to the lower atmospheric layers, v/here the abimdance of H 2 O is 
relatively high. If this oxidation mechanism is correct, then 
the turbopause on ^krs must be at least as high as the region of 
photochemical activity, (i,e, at levels above 75 km), 

V/hile it is fairly certain that the photochemical equilib- 
ri\im distribution of gases shown in Figure 15 is not a correct 
representation of the actual distribution of gases in the 
Martian atmosphere, no means is available at present for making 



a more accxirate determination of the chemical structure, due to 
uncertainties with regard to the principal recombination reac- 
tions involved, and the relative importance of eddy and molecu- 
lar diffusion in the region of photochemical activity. 

It is possible, hov/ever, to set reasonable upper and lower 
limits to the total number density at higher levels, based on 
the photochemical equilibrium distribution, for each model 
considered. The effect of strong eddy mixing in the upper 
layers would be an upward transport of CO 2 compensated by a 
dovmward flux of CO and 0. Such a rearrangement of the constit- 
uents would imply that the total number density near the top of 
the models in the presence of strong mixing would be the same as 
indicated by the 0 - CO curve (Figure 15) which, under photo- 
chemical equilibrium conditions, falls off with the scale height 
of C 02 « 'hhe other hand, molecular diffusion under the influ- 
ence of Martian gravity would tend to increase the concentra- 
tions of all constituents at higher levels over the values 
indicated in Figure 15, with each gas being distriouted according 
to its own scale height above the turbopause. Since the concen- 
trations of the lighter dissociation products 0 and CO would 
decrease more slowly with altitude than is indicated by photo- 
chemical considerations, one would expect these gases to 
rapidly become more abundant than the 0 - CO profile in Figure 
15 indicates, A reasonable upper limit to the total number 
density at levels some distance above the region of maximum 
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photochemical activity would be given by the diffusive profile 
for 0, 

The total number density determined from basic thermal 
considerations alone, neglecting photochemistry, mixing and 
molecular diffusion may, therefore, be thought of as a lower 
limit to the actual number density at higher levels, while the 
diffusive 0 profile above the level of maximum 0 concentration 
may be considered an upper limit* An estimated lower limit to 
the neutral particle number density at ionospheric levels in 
Models I thru VII was obtained as follows. The temperature at 
the top of each of these thermal models was adopted as character- 
istic of the mean temperature between the 100 km level and the 
level of the ionosphere (124 km). Using this temperature, an 
estimate of the mean scale height for CO 2 ( or a CO 2 - mix- 
ture), in this region of the atmosphere was obtained for each 
model, and the 100 km number density was extrapolated to the 
ionospheric level. The lower limit to the nuraber density just 
below ionospheric levels (at 120 km) has already been cited 
for Models VIII and IX, so that a similar extrapolation yields 
the number densities for these models at the ion peak. The 
upper limit to the neutral particle density in the ionosphere 
was obtained by extrapolating the 0 concentration from the 
level of the 0 maximum to the 124 km level, using a character- 
istic scale height for 0 determined from the basic thermal struc- 
ture for each model. The upper and lower limits to the total 
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neutral number density at 124 km obtained in this v/ay are 
listed in Table V. 

If the computed thermal structures for the nine basic 
models (Figures 11, 12 and 13) are accepted, then the lov;er 
limits given in Table V for the neutral particle density at the 
level of the ionosphere in each model is indeed a lower limit. 

No physical process can produce a number density smaller than 
this value for the given thermal structure. The thermal pro- 
files depicted in the figures may, however, constitute slight 
overestimates (by /v 5®K) of the true radiative equilibrium 
temperatures in the upper atmosphere for each case, due to over- 
emphasis of both solar energy input and absorption of weak- 
interval surface radiation in the higher layers of the models. 

If, in fact, solar radiation has negligible influence on the 
atmospheric temperature structure belov; ionospheric levels, and 
if the temperature maximum indicated in the figures is not a 
real feature of the Niartian atmosphere , then both the lower and 
upper limits to the neutral ionospheric number density given in 
Table V v/ould be reduced by a/ 50>j, 

If one accepts a surface temperature in the range 1S0"K to 
200*K, a surface pressure in the range 5 nib to S mo and a compo- 
sition of to lOO'p CO 2 as an accurate description of the 
basic conditions in the lower llartian atmosphere at the time and 
place of lariner IV occultation, then these considerations indi- 
cate that the neutral number density in the observed ionosphere 
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TABLE V 

ESTIKATED RANGE OF NEUTRAL NUI-BER DENSITY AT ICNCSPHSRIC LEVEL 


Model 

Lower Limit 

Upper Limit 


(raol/cm^) 

(mol/cm^) 

I 

3.6x10^ 

1.3xl0ll 

II 

1.6x10^® 

2.3xl0ll 

III 

3.9x10^° 

3./pcloll 

IV 

1.2x10^^ 

5.0xloll 

V 

7.6x10^ 

1.7xloll 

VI 

2,3x10^° 

2 . 5x10 1^ 

VII 

7.2x10^° 

4.2xloll 

Till 

1.6x101® 

2.Axloll 

IX 

3.7xl0l® 

2.9xl0ll 
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10 3 11 3 

was betv;eon a/ 10 mol/cm and ^ 10 mol /cm « The smaller 
value is representative of mixed conditions, with CO^ the most 
abundant gas at ionospheric levels, and the larger value is 
characteristic of a graviationally separated upper atrao sphere, 
v:ith 0 the principal constituent in the ionosphere, nil of the 

9 / 3 

values within this range are much larger than the 2 x 10 mol/cra 
proposed by Johnson (1965) and more than an order of magnitude 
smaller than the 5 x 10 mol/cm*^ obtained by Chamberlain and 
1-icElroy (i960) , using the lower atmospheric model of Prabhakara 
and hogan (1965) • On these grounds, the ionospheric models 
proposed by these authors cannot be correct. 

Ionosphere , 

In view of the results of tl.e photochemical and radia.tive 

transfer calculations described above, it appears that ti.e most 

abundant gas at ionospheric levels on ^lars must be either CC2 

or 0, the former if mixing governs the distribution of gases and 

the latter if the constituents are gravitationally separated 

at that level. Moreover, indications are that tl,e total number 

density at the height of t'ne ionosphere, if the turbopause is 

1 0 3 

located near ionospheric levels v;ould be a/ 10 mol/cm , while 
if the turbopause is found nearer the region of maximum photo- 
chemical activity, the total number density in the ionosphere 
v;ould be A/ 10^^ mol/cm^« In the present section, an attempt 
is made to determine the principal constituent, main ionic 
reactions and atmospheric conditions at the level of the 
ionosphere on Mars, 
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Suppose that 0 is the main constituent at and above the 
ionospheric level on fers, and is distributed with its ow-n scale 
height. Then 0"^ will be the principal ion formed by the ionizing 
EUV solar radiation, v;ith a maximum number density of 0"^ occurr- 
ing at the level where the neutral 0 density n(0) is given by 


cos© 

n(0) = 

a(0) H(0) 


(23) 


where 0 is the solar zenith angle, 0(0) is the mean photoion- 
ization cross section for 0, and H(0) is the scale height for 0 
at the ion peak. The value of 0 corresponding to I-iiriner IV 
immersion is a/ 70® , the mean photoionization cross section for 0 
is /V 10’ cm^ (Hinteregger, ^ al., 1965), and if the thermal 
models calculated in the present study are accurate, the scale 
height of 0 at ionospheric levels is /v 20 km. Thus we obtain 
a value of 2 x 10^® mol/cm^ for the 0 density at the level 
where a maximum in O’*” production occurs. The rate of formation 
of o’*" ions at the maximum is given by 

Q(0'^) =FCr(0)n(C) exp(-l) (24) 

max 

v/here F is the total ionizing flux at the top of the atmosphere, 

10 2 

'v 10 photons/cm'^sec (Hinteregger, at a^. , 1965). The value 
of Q(0 ) at the ion peak is therefore 7.4 x 10^ ion pairs/cm^se c, 
and if a balance between local production and loss of ions is 
assumed, 

Q(0*) - o(n2(e) « 0 (25) 


where is an effective recombination coefficient for electrons 
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with ions. Since a maximum ion number density of /v 9 x 10 
electrons/cm"^ v;as observed by tferiner IV, an estimate of about 
9.1 X 10“ cmi sec for oi is obtained from this balance relation- 
ship. Such a value is characteristic of the dissociative recom- 
bination coefficient for the reaction 


© 2 ^ + e — » 0+0 


(26) 


(Biondi, 1964), and suggests that a faster reaction rapidly 
transforms the EUV produced 0^ into 02^ » Since it is well known 
that an ab\jndance of CO 2 is present in the Martian atmosphere, a 
likely mechanism for this transformation is 

O"^ + CO2 —*00+02'*’ , ( 27 ) 

for which a rate coefficient of 1.2 x 10 ^ cm^sec ^ has been 
proposed (Feshenfeld, et aj., , 1966). 

At higher levels, where the atmosphere is optically thin 
in the ionizing frequencies, the ion production rate is given by 

Q(0'^) = F (T(0)n{0) , (2S) 

or, from the above values of F and C(0), 

Q(0'^) = 10“*^ n(0) . (29) 


On the other hand, the loss rate of 0"^ from reaction (27) is 


L(0'^) = 1.2 X 10"^ n(0'*‘) n(C02) 


(30) 


A balance between production and loss of 0 ions implies that 
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lO"*^ n(0) - 1,2 X 10"^ n(0'*") nCCOg) 


-9 


so that 


n(O^) S’ 10^ 


n(0) 

n(C02) 


(31) 

(32) 


The suppression of a diffusion controlled F2 ion peak to a 
non-detectable level (n(0 ) less than 10^ mol/cm^) requires 
that the abundance of C(^ remain large enough to rapidly remove 
0 ions via reaction (27). From expression (32), if n(O^) is to 
remain less than /w 10^, the inequality 


n(0) 
n(C02) 

must hold, up to the level at which diffusion becomes the 

dominant mechanism for 0^ removal. 

At levels in the atmosphere where local loss of 0^ ions is 

governed by ambipolar diffusion, the loss rate L(O^) may be 

written . 

+ + D(0 ,0) 

L(0"^) = n(0"^) ^ , (34) 

H(0)2 

where D(0"'^,0) is the ambipolar diffusion coefficient for 0^ ions 
in atomic oxygen, inversely proportional to n(0). The level at 
which diffusive loss of 0'*’ is equal to the loss via reaction (27) 
is reached at the point where 

Q 0(0”^, 0) 

1.2 X 10“9 n(C02) - . (35) 

^ H(0)2 

Donahue (1966) finds that if H(0) is ^50 km at this level 


< 10 " 


(33) 
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(T ^ 360® K), then the CO^ concentration is related to the 
n(0) concentration via equation (35) hy 

n(CO^) - 3 X 10^ / n(0) , (36) 

at the level where ambipolar diffusion becomes important* 

Since inequality (33) must hold at this level also, 

simultaneous solution of (33) and (36) results in a value for 

A *1 ^ 

n{C0^) of 2 X 10^ and a value for n(0) of less than a/ 2 x 10 

mol/cm^ at the diffusion level. Such a value for 0 requires 
a decrease of two orders of magnitude (4*5 scale heights) in 
n(0) between the ion peak and the diffusion level. Due to the 
smaller scale height of CO^ (2,75 CO 2 scale heights per 0 scale 
height), a decrease of n(C 02 ) by more than six orders of magni- 
tude (12*4 scale heights) between ion maximum and diffusion 
level is implied. Thus a number density of CO 2 of ^ 2 x 10 
mol/cm^ at the diffusion level demands a number density of 
/V 5 X 10^^ at the ion peak. Such a value is larger by an order 
of magnitude than the number of 0 atoms initially assumed to be 
present as the main neutral constituent near the ion peak. 

Thus, an irreconcilible contradiction arises from the assumption 
of a predominately 0 - — 02^ ionosphere* The concentration 

of CO^ required at higher levels to suppress an F2 peak, implies 
that CO 2 must be present near the ion peak in an abmdance 
greater by far than that of 0 required to produce the peak. 
Suppose, on the other hand, that CO 2 is the dominant 
constituent at 124 km, but that at higher levels diffusive 


117 


separation occurs and 0 rapidly predominates. Then the main ion 
produced at 124 km would be CO^ , The mean ionization cross- 
section for CO^ is slightly larger ( 2x) than the cross-secticn 

for 0 (Hinteregger, ^al., 1965), but the smaller scale height 
(by a factor of ^2,75) implies (from the CO^ analogy of equaticn 
(23)) that the maximum production of ions would occur where 

the neutral density of CO^ is 4. 2.3 x 10^^ raol/cm^. The maximum 
rate of CO^ ion production would be (from the C0« analogy to 
equation (24)) a/ 1,7 x 10 ion pairs/cm^ sec, A balance between 
production and loss of CO^ ions (via the CO^ analogy of equation 
(25)) results in an estimate of the effective recombination coef- 
ficient for ^^2 ions of ^ 2,1 x 10 cm^ sec”^. Such a value 
is slightly smaller than the 3.S i 0.5 x lO"*^ cm^/sec estimated 
by heller and Biondi (1967) for dissociative recombination of 
^^2 with electrons. Nevertheless, it suggests that this 

may be the chief mechanism for ion removal. It is only reason- 
able to expect that, even if 00^ is the main constituent in the 
ionosphere, at higher levels 0 will eventually predominate and 
0 will be the dominant ion formed. The same conditions must 
therefore be met for the suppression of the F2 ion peak as 
described previously, leading to the conclusion that at the 
level where the CO^ number density is a/ 2 x 10^, the number 

g 

density of 0 atoms must be less than a/ 2 x 10 . In this case, 
a decrease by a factor of 10^ (9 scale heights) in CO^ concen- 
tration must occur between ion peak and diffusion level, 
corresponding to a decrease in the 0 concentration by a factor 
of 2»6 X 10^ (3*27 scale heights). Thus, the 0 number density 
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at the ion peak would be at most 4,5 x 10 mol/cm , well belov; 
the CO2 number density initially assumed to be present at that 
level, so that, with this interpretation, no contradiction arises, 
A simple exercise demonstrates the ability of this inter- 
pretation to produce the ionospheric profile observed (Figure S), 

It was assumed that the concentration of CO at the ion peak 

10 3 ^ 

was 2 X 10 mol/cm , that the 0 concentration at this level is 

g 3 

lOx smaller (2 x 10 mol/cm ) , and that both constituents are 

distributed according to their own scale heights above the peak. 

It was further assumed that both constituents are ionized by a 

10 2 

solar flux of 10 photons/cm sec, CO^ with a cross-section of 

2 X 10“^*^ cm^ and 0 with a cross-section of 10”^^ cm^, Disso- 

+ + 

ciative recombination was assumed to proceed for CO^ and 0^ 
with the rates of Weller and Biondi (1967) and Biondi (1964)» 
CO^^, 0^^ and total ion number densities were calculated at unit 
scale heights above the peak. The resulting neutral particle 
and ion densities are listed in Table VI, It can be seen that 
the total ion number density at the peak is quite close to the 
observed value, with approximately two-thirds of the ions at the 
peak CO^^ and one-third 0^^; at higher levels, the ratio O2 
increases, and at a level six CO2 scale heights above the peak, 
this ratio is -v 3* By relating the CO^ scale height intervals 
to altitude above the ion peak ( 124 km) , it is possible to fit 
these total ion number densities to the observed electron 
profile (Figure d), A reasonable fit is obtained using the 
tabulated values, if a scale height gradient of 1 km/scale 


TABLE VI 


IONOSPHERIC MODEL 


Height 

Above 

Ion 

Peak 

CM 

0 

0 

Number Density 
C02^ 0 

(cm~^) 

V 

Electrons 

0 

2.0(10)** 

6.5(4) 

2.0(9) 

2.7(4) 

9.2(4) 

1 

7.4(9) 

5.3(4) 

1.4(9) 

3.0(4) 

g.3(4) 

2 

2.7(9) 

3.6(4) 

9.7(g) 

2.g(4) 

6.4(4) 

3 

9M&) 

2.4(4) 

6.7(g) 

2.6(4) 

5.0(4) 

4 

3.6(g) 

1.4(4) 

4.7(g) 

2.2(4) 

3.6(4) 

5 

1.4(g) 

0.9(4) 

3.2(g) 

l.g(4) 

2.7(4) 

6 

5.0(7) 

0.6(4) 

2.2(g) 

1.4(4) 

2.0(4) 


Height* 

Above 

Ion 

Peak 

Estimated 

Altitude 

(km) 

Estimated 

Temperature 

CK) 



0 

124 


140 



1 

131 


160 



2 

139 


igo 



3 

14g 


200 



4 

153 


220 



5 

169 


240 



6 

Igl 


260 



measured In scale heights of CO2 

powers of 10 are enclosed in parentheses, e.g., 

2.0(10) - 2.0 X 10-^. 
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height CO^ is used above the ion peak. This scale height 
gradient is consistent with the thermospheric model of 
Chamberlain and McElroy (1966), The altitute and temperature 
corresponding to each ion density level for this fit are also 
given in Table VI, 
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SUMMARY 

The present construction of model atmospheres for Mars 
from detailed considerations of radiative transfer demonstrates 
that the low atmospheric temperatures proposed in order to 
identify the observed ionosphere as an F2 layer (Johnson, 1965j 
Fjeldbo, ^ al*» 1966, a, b) cannot be attained if the physics 
of the atmosphere is properly taken into account. The present 
calculations also indicate that, if values of the basic atmos- 
pheric parameters more representative of the time and place of 
Mariner IV occulta tion are adopted, and if vibrational relaxa- 
tion of COg and Doppler broadening of the spectral lines are 
considered, then the Martian atmosphere is probably significantly 
colder than was shown by Model I of Prabhakara and Hogan (1965)# 
Thus, the E- layer interpretation of Chamberlain and McElroy 
(1966) which was based upon this earlier model of the lower 
Martian atmosphere has no solid foundation. 

For lower atmospheric conditions suggested by theory and 
observations, it is found, from considerations of photochemical, 
mixing and molecular diffusion processes in the light of the 
thermal structures obtained, that probable lower and upper limits 
to the neutral number density in the observed Martian ionosphere 
are ^ 10^® mol/cm^ and /v 10^^ mol/cm^ respectively, the former, 
CO 2 , being attained if mixing extaids to near-ionospheric levels, 
and the latter, 0, being reached if gravitational separation 
predominates above the photochemically active region. 
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For the calculated thermal structures, at a zenith angle of 

+ + 

'V 70”, the maximum rate of production of either CO^ or 0 ions 

occxurs at the level where the neutral particle number density is 

A' 2 X 10^® mol/cm"^. The effective recombination rate required 

to produce the observed electron number density maximum is 
—7 3 

'V 10 cm^/sec, characteristic of dissociative recombination of 

+ + + + 

O2 t or possibly, CO^ • The interpretation of an 0 - 0 " 

ionosphere requires that CO^ be the main constituent near the 

ion peak, in order to account for the suppression of a diffusion- 

controlled ion peak at higher levels, and thus involves a 

contradiction* On the other hand, the assumption that CO^ and 

CO ^ are the main neutral and ionized constituents near the ion 
2 

peak results in no apparent contradiction. In this case, the 
criterion for suppression of an F2 peak at higher levels requires 
that the 0 concentration near the ion peak be at least an order 
of magnitude smaller than the local CO2 concentration. 

Calculations from simple ionospheric theory of the CO^ and 

0 abundoOices at and near the ion peak at 124 km, assuming 

2 1 0 3 Q 3 

neutral densities of 2 x 10 mol/cm CO^ and 2 x 10 mol/cm 0 

at the peak, show that, although the concentration of 0 is only 

one-tenth of that of CO2 at 124 km, nevertheless O2 ions 

constitute about one-third of the total ion population at that 

level, due to the slower recombination rate for these ions* At 

a level three CO2 scale heights ( /v 25 km) above the peak, CO2 

and 02^ ions are equally plentifxxl, while at higher levels where 

the CO2 concentration falls off, 02'^, and eventually 0*, 
becomes the dominant ion* 
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The probable structure of the Tertian atmosphere and ionos- 
phere at the time and place of occultation indicated by the 
present investigation may be described as follows. The surface 
temperature was within the range 160®K to 200* K, the surface 
pressure was between 5 mb and S mb, and the composition of the 
atmosphere was dO^ to 100^ ^^2* tropopause level lay between 

2 km and 3# 5 km above the surface of I4ars, with temperatures 
decreasing at higher levels to values of 140* K - 150° K, 13 5° K - 
U5®K, and 130*K - 140°K at 50 km, 100 km, and 125 km respec- 
tively. The atmosphere is mixed up to near-ionospheric levels, 
but at and above the obs^ved ion maximum molecular diffusion 
under the influence of Martian gravity is likely to predominate. 
There is little dissociation of CO^ below ionospheric levels. 

The observed ion peak on I^hrs is an FI- layer, with dissociative 
recombination of molecular ions the chief mechanism for electron 
removal* The neutral particle density at the ion peak, nainly 
due to COgf is 2 X 10 mol/cm*^. Two-thirds of the ion popula- 
tion at the peak is C0«^, one-third is 0 A few CO, scale 

^ ^ ^ 

heights above the peak, Og becomes the dominant ion, and at 
still greater altitudes, 0^ ions must predominate. A gradient 
in scale height (temperature) of about 1 kiVC 02 scale height 
exists above the ion peak. It is possible, and, indeed, likely, 
that a true thermospheric structure, such as found by Chamberlain 
and KcElroy (1966) and Stewart (1967), is present in the upper 
atmosphere of Mars* 
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The above description is consistent with the I'ariner IV 
results as well as with recent spectroscopic evidence. In 
addition, it is in agreement with careful studies which 
indicate the presence of a warm thermosphere on ^^ars. 
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APPENDIX A 

THE EQUATION OF RADIATIVE TRANSFER 


The basic equation describing the transfer of radiation 
through a plane parallel atmosphere, Schwarzschild*s equation, 
may be vrritten 

A — - j(v,t*) (A.i) 

where I(v,yi|,1: ) is the specific intensity of the radiation 
is the soxirce function 
V is the frequency of the radiation 
/W “ cos © 

© is the zenith angle 

't* is the optical depth (d't* = - n(z)k(v,z)dz) 
k(v,z) is the absorption coefficient 
n(z) is the number density of absorbing molecules 
z is the geometric height above the planetary surface. 

If suitable boundary conditions are specified, the equation of 
radiative transfer (A.I) can be solved by standard methods for 
the upward and downward directed specific intensities at optical 
depth t-”^, I^(p,^,'t*) and respectively. Physically 

meaningful boundary conditions for a planetary atmosphere require 
that the downward intensity at the top of the atmosphere (t = 0) 
be zero, 

l^(v,^,0) = 0 , (A. 2) 

and that the upward directed intensity at the lower boundary 

) be the blackbody intensity of the planetary surface, 

6 
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(A. 3) 


One method of solution of equation (A.l), involving Green’s 
functions, proceeds as follows. 

Equation (A.l), an inhomogeneous, linear, first order diffe- 
rential equation, has a solution since the corresponding homoge- 
neous equation, 


h 


dl(v,//,qr*) 

d't* 




(A.4) 


has a solution. The solution of (A.4) for the upward intensity 
is of the general form 

* C* exp(^t*^) (A. 5) 

where C’ is an arbitrary constant. Then the solution of (A.l) 
may be expressed as 




H('t*,t)r(t)dt 


(A.6) 


where C” is a constant to be determined by the lower boundary 
condition, and4>(t'*) is bhe solution of (A.4), 

q)(tr‘‘) » C’ exp(tr*/^) • (A. 7) 

Also, r(t) is the inhomogeneous terra in the transfer equation, 

r(t) = - J(v,t)//< , (A. 8) 

and H(t*,t) is the one-sided Green’s function, here given by 


H(^*,t) « -4>(t;*)/vy (t) - - exp(-(t-^*)//<) . (A. 9) 

Substituting equations (A. 7), (A. 8) and (A. 9) into (A.6), we 
obtain 

I^( Vi^iTr") - C’C”exp('t*y^) + /* ® J(V,t)exp^(t-'b 

^ (A. 10) 
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Applying the lower boundary condition (A. 3) we have 


so that 


- 0-C<<exp(t*g//,) - B(v,t* ) , (A.ll) 


C>0" . B(v,-t:‘g)exp(-t*g/y<) 


(A. 12) 


Substituting from expression (A. 12) for C’C»» in (A. 10), the 
relation for the upward directed specific intensity becomes 


/r* / . \ (A. 13) 

^ ^ J(v,t)exp^(t-'^'‘)/^^ . 

By a similar sequence of operations, we derive the expres- 
sion for the downward directed specific intensity. The solution 
of (A, 14) for (v,y/,*tr*) is of the same form as given by (A.5)« 
Then the general solution of (A.l) will be of the type 

I*(v,//,t*) =C»»^;('t*) + H(t^t)r(t)dt (A. 14) 

where C * ’ is a constant to be determined by the upper boundary 
condition, v|J ('t*) is given by (A. 7), and r(t) is given by (A,S), 
H(t^',t) is here given by (A. 9). Substituting from (A. 7), (A.S) 
and (A. 9) into (A.14)> we obtain 

» C*C*»exp(tr*/^) + J(v,t)exp(^(t-'l?‘)^^ 

At the upper boundary, (A, 2) applies, so that 


I^(V,yV,0) = C»C»» « 0 


(A. 16) 


Hence, expressicm (A. 15) for the downward directed specific 
intensity becomes 



0 


\ V H 

For 9 between - and /f ^ positive value, so that 

upward directed intensities will be regarded as positive. For 0 

between"^ and M is negative, and the downward directed 

2 2 " 

intensities will therefore be negative. 

The upward and downward specific intensities may be v/ritten 
in terms of the atmospheric transmission function 't'. The trans- 
mission of the atmosphere at frequency V and zenith angle arc cos 
between optical depths and t is 

Hence, 






(A. IB) 
(A.19) 


and equations (A. 13) and (A. 17) may be written 

**) « ^) 'CtV,/f,'Cg-ir J 

* 


-U — j; — “ 


(A. 20) 


and 


dt 


(A. 21) 


J(v.t) 

respectively* 

For the present problem, the specific intensities were 
expressed in terns of 0 and z, rather than in terms of ^ and ^ 
The relationship between tr and z is given by 


n($)k(V,5)d5 


(A.22) 
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where z.p is the top of the atmosphere. (The actual height 
chosen for the upper boundary is arbitrary, provided no signifi- 
cant contribution to the net flux of infrared energy originates 
above this level.) The maximum value of 'Ir in the atmosphere 
occurs at the lower boundary z 

o 

•r*g(2g) - n(s)k(y,s)ds . (A. 23) 

s 

V/ith these considerations, the specific intensities (A. 20) and 
(A .21) may be written 


I^V,e,z) 

and 


B(V,Zg) t’(v>,e,z-Zg) 



J(v,s) 


^'^(v,e,z-s) 


(A.24) 


:^(v,e,z) 


J(v,s) 


SI— 


dS 


(A. 25) 


corresponding to solutions of the transfer equation (A.l) in 
the form 
dl (v ,©,z) 


ds 


= - n(z)k(v>,z)I(v,e,z) + n(z)k(y,z) J(V,z) 


(A. 26) 


where ds = sec dz 
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APPENDIX B 

VIBRATIONAL RELAXATION OF THE CARBON DIOXIDE MOLECULE 

An individual molecule may possess four modes of internal 
energy: translational, vibrational, rotational and electronic. 

Only the first three of these energy modes will be considered 
here. Because of the relatively large energies involved in 
electronic transitions, thermal radiation from excited electronic 
levels is not, in general, important in planetary atmospheres. 

The translational or kinetic mode of molecular energy is 
unquantized; a molecule may possess any amount of kinetic energy 
and may freely exchange this energy in collisions with other 
molecules. (A molecular collision, by definition, involves an 
interchange of kinetic energy.) Since there is a rapid adjust- 
ment to equilibrium (Boltzmann's law) in this mode, the molecular 
velocity distribution is continuous and Maxwellian in nature up 
to very great heights in a planetary atmosphere (Spitzer, 1952). 
The mean translational energj'’ of the molecules is an observable 
quantity and therefore the temperature of a gas is usually 
defined as the kinetic temperature or that associated with the 
molecular velocities. 

The vibrational and rotational internal energy modes are 
quantized; only discrete energy levels are possible and a number 
of rotational energy levels are associated with each vibrational 
level. Thus, a given vibrational transition can occur in 
combination with any one of a number of rotational transitions. 
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The probability of energy transfer between vibration - rotation 
and translation during a collision is less than one, and a 
number of collisions f normally occurs before an energy exchange 
between vibrational - rotational and translational modes takes 
place. If the collision rate in the gas is f^, then the relaxa- 
tion time \ , or mean time for loss or gain of energy by the 
vibration - rotation modes by collision is given by 

f 

'^•-7 (B.l) 

where \ is a function of the pressure p and the temperature. 

Since the collision frequency is directly proportional to the 
pressxire, at a given temperature 


(B.2) 

where is the relaxation time at a reference pressure p . The 
pressure variation of X is the one of importance, since pressure 
varies by many orders of magnitude in a planetary atmosphere, 


while temperatures are ordinarily confined to a relatively 
narrow range. 

In addition to excitation and deexcitation by collisional 
processes, transitions between vibrational or rotational levels 
may result from the absorption or emission of a quantum of 
radiant energy. Any excited molecule has a (statistical) natural 
lifetime for spontaneous deexcitation by a radiative transition, 
Curtis and Goody (1956) have shown that, if the Planck function 
does not vary much within a vibration - rotation band, a 
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radiative lifetime for such a band may be expressed as 


9 


■1 - 


/k^dv 
\ c / y ^ 


(B.3) 


where Cp is the radiative lifetime, Vq is the central frequency 
of the band, c is the velocity of light and ^kyd^ is the band 
intensity or integrated absorption. 

V/hen pressure is high and collisions are frequent there is 
a rapid exchange of energy between vibration - rotation and trans- 
lation, and the energy levels of the vibrational and rotational 
states are populated in a Boltzmann distribution due to colli- 
sions, Under these circumstances, the vibrational or rotational 
temperature is equal to the kinetic temperature. The gas is 
said to be in a state of local thermodynamic equilibrium (LTE) 
and it radiates according to Kirchhoff’s and Planck’s law. Due 
to the frequent collisions, there is strong coupling between 
matter and the radiation field. 

As pressure decreases, collisions become less frequent and 


the transfer of energy between vibration - rotation and trans- 
lation becomes less rapid. The relaxation time increases and 
eventually becomes longer than the radiative lifetime, V/hen this 
occurs, the energy involved in deexcitations will be lost to the 
radiation field. Moreover, when a photon with a frequency lying 
within the range of the relaxing band is absorbed by a molecule, 
it will on the average be reemitted without passing into the 
kinetic energy of the gas or translational mode. Under these 
conditions, the gas is said to undergo ’’relaxation" and the 
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atmosphere is said to be "scattering”. Collisions are insuffi- 
cient to maintain the vibration - rotation energy levels in a 
Boltzmann distribution in the presence of the radiation field. 
The radiation field will then determine the populations of these 
energy levels, and the coupling between matter and radiation 
will be weak. The gas no longer obeys Kirchhoff*s or Planck's 
laws and is said to be in a state of non-local thermodynamic 
equilibrium (non-LTE) , Therefore, at low pressures, one cannot 
properly relate the temperature of the gas to the vibrational - 
rotational energy levels without considering the relaxation 
process. 

The 15^ vibration - rotation band of is certainly the 
most important band for the transfer of infrared radiation in 
the lower atmosphere of Mars, Vibrational relaxation in this 
band must be taken into account in any radiative transfer study 
since it occurs at altitudes where CO^ is present in abundant 
quantities. The level of vibrational relaxation in the I5 ^ COg 
band depends simply on the adopted values of the relaxation time 
and radiative lifetime. 

Rotational relaxation of CO^ will be of no importance in 
this study. The energies involved in transitions between rota- 
tional energy levels is comparatively small, so that the proba- 
of a rotation - translation transfer of energy in 
collisions is much greater than the probability of a vibration - 
translation interchange, xne rotational relaxation time will 
therefore be many orders of magnitude smaller than the 
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vibrational relaxation time (Lambert, 1962), and rotational 
relaxation could occur only at very high altitudes (low pressures) 
where COg is likely to be dissociated. 

The 15 ^ vibration - rotation band of CO 2 is composed of a 
fundamental (v >2 or ’’bending” mode of vibration) at 667«40 cm 
and thirteen overtone and combination bands. The fundamental 
is by far the most intense of these bands, accounting for about 
90 percent of the overall 15^ band intensity (Yamamoto and 
Sasamori, 195S; Madden, 1961). The band intensities of the sub- 
bands in the 15 region are much smaller (Young, I 964 ). Accord- 
ing to equation (B.3), therefore, the radiative lifetimes of the 
sub-bands are considerably longer than the lifetime of the 
fundamental* Unfortunately, no estimates of the relaxation 
times of the CO 2 sub-bands near Vj>jA have appeared in the liter- 
ature, If however, the relaxation time of these bands is of 
the same order as the relaxation time of the fundamental (Young, 
1964 ), they will relax at much higher levels, due to their 
longer radiative lifetimes. One might then assume that these 
subsidiary bands are in LTE in the lower region of the i^artian 
atmosphere. 

Since the contributions of the overtone and combination 
bands to the infrared energy flux are proportional to their 
intensities, their influence on the infrared transfer is small 
compared to the influence of the fundamental. Therefore, it is 
assumed in the present study, that the sub-bands relax in the 
same manner as the fundamental. The radiative lifetime and 
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relaxation time for the 15^ fundamental is used throughout the 
15yf region. This assumption permits an enormous simplification 
of the computational procedure involved in calculating a temper- 
ature structure from energy balance considerations. 

According to Madden (1961) the band intensity of the 15 
CO^ fundamental ts m cm’^ atm*!. Using this value for the 
intensity, one obtains a radiative lifetime of 4J2 x 10**^ sec. 
Experimental determinations of the vibrational relaxation time 
for the fundamental (Herzfeld and Litovitz, 1959; Lambert, 1962; 
Marriott, 1966) give estimates ranging from 10”^ sec to lO”^ sec 
at a pressure of one atmosphere. Vibrational relaxation in the 
15 jA CO 2 band should then occur at an altitude where the pressure 
is between /v 2,5 x lO”^ mb and 2,5 x 10“^ mb. 

The \ - 6^ near infrared CO^ bands also undergo vibration- 
al relaxation at low pressures. These bands are important for 
the absorption of solar energy, and this absorption could become 
large at high levels where the lines are unsaturated. However, 
a study of solar heating in the Earth »s atmosphere by Me Clatchey 
( 1966 ) has shown that vibrational relaxation in these bands 
probably occurs far below the levels at which the stronger lines 
become unsaturated, so that the solar heating rates due to these 
bands are probably very small throughout the Earth’s atmosphere. 
The same arguments would apply to Mars, and solar near infrared 
heating in the Martian atmosphere must be small also. According 
to Me Glatchey (1966), relaxation of the intense V.* band of CO 

2 

near 4,3^ begins to exert a strong damping effect on the 
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heating in this band at pressures between 10 mb and 10 mb, 
and at lower pressures, almost complete scattering of solar 
energy occurs. 

Thus, the near infrared bands relax at higher pressures 
than does the 15 00^ band, and their influence on the tempera 
ture structure should be of little consequence. In order to 
assess their effect on the equilibrium distribution of tempera- 
ture, however, a mean relaxation time X' and radiative lifetime 
(p' may be assigned to these bands as parameters for variation. 
The fraction of the interacting solar energy in these bands 
which is not scattered, but which passes into the kinetic mode 
and appears as local atmospheric heating will then be given by 
( 1 ^ where, as in the case of the 15^ band, the 

relaxation time is inversely proportional to pressure. 
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APPENDIX C 

RADIATIVE TRANSFER UNDER LTE AND NON-LTE CCNDITICNS 


The integral form of the transfer equation (A. 26) 


dvdu> 


- JJ »e»z)n 


(z)k(v,z) dvdu> 


(C.l) 


where <o is solid angle (du> «= 2trsined0), relates the divergence 
of the net flux of infrared radiation to the difference between 
the absorbed and emitted energies per unit volume at any level z* 
Here, frequency integration is performed over the transferring 
band(s) and integration with respect to solid angle cc is perform- 
ed over both hemispheres. 

If local dynamical or chemical sources and sinks of energy 
are absent, and if no radiation is absorbed at frequencies beyond 
the range of the integration in (C.l), radiative equilibrium 
conditions require that the infrared flux divergence be zero, so 
that there is a balance between absorbed and emitted infrared 
energies 



I(V,©,z)n(z)k(v,z)d>>dco 



J(V ,z)n(z)k(v ,z)dwdo 


0 


(C.2) 


Since I(v,©,z) is itself a function of J(v,z) ((A. 24 ) and (A. 25)), 
expression (A,26) is an integro-differential equation for J(v,z), 
which may be used to determine an equilibrium atmospheric 
structure, provided a relationship can be established between 
J(v,z) and the local atmospheric conditions (including the 
local radiative field). In order to obtain a unique thermal 



structure, however, the surface temperature T (which determines 
B(v,z_)) and a form for the transmission fiinction ^(v ,6 , jz- Sj ) 

o 

must be specified. The form and method of evaltiation of the 
infrared transmission function is described in Appendix D. 

If sources and sinks of energy outside of the range of the 
frequency integration in (C.l) are present, then, under radiative 
equilibrium conditions, the local infrared flux divergence must 
be equal to the local energy gain or loss Q(z) from these sources 
and sinks: 


V/^e then have 



dl(v,e,z) 

ds 


dvdo) 


Q(z) 


(C.3) 


Q(z) + 



I(V,©,z)n(z)k(V,z) d-Jdoi 

J(V ,z)n(z)k(v |Z) dvd<o “ 0 


-//' 


(C.4) 


which, with the above considerations, may be solved for a unique 
atmospheric structure also, if Q(z) is known. 

The relationship between J(v,z) and the local conditions in 
the atmosphere varies with the collision frequency. In regions 


where collisions are frequent, J(->>,z) will simply be the ’’Planck” 
emission, expressed in terms of the gas temperature, and given by 
the relationship • 


J(V,z) = B(V,T(z)) 


(C.5) 


2hv^ 

c^(exp(hv/KT(z) ) - 1) 
where h is Planck*s constant, K is Boltzmann’s constant and c is 
the velocity of light. As discussed in Appendix B, in low 
pressure regions where ’’relaxation” processes occur, molecular 
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collisions are too infrequent to mintain a Boltzmann distributicn 

among the vibration - rotation energy levels. The atmosphere in 

such a region does not. therefore, radiate according to (C.5) and 

another form is required to properly relate J(v,z) to local atmos- 
pheric conditions. 


In their investigation of the problem of radiative transfer 
in a vibrationally relaxing gas, Curtis and Goody (1956) derived 
a general expression for the emission of the gas J(v,s) which is 
applicable under low pressures. They have shown that if Av , the 
band width is narrow, so that variation of the Planck function 
B(v,z) over the band is small, then J(v,z) will have the same 
dependence on frequency as B(V,s) and will be isotropic. From 


their equation (10) (in the present notation), 

J(z) = B(z) + ^ 

9 n(z)k(v,z)dV 

where J(z) and B(z) are the values of J(v,z) and B(v,z) at the 

band center (J(V,z)a J(z) and B(V,z) B(z) ) , and where 

dl(v,6,z) 


(C.6) 


h(z) 


-// 


ds 


dvdco 


it follows that the source function is related to the Planck 
function by 


. /' \ ff d^doo \ 

J(v ,z) = B(v,z) f 1 - — 

\ 9 yy^(V,z)n(z)k(v,z)dvdc«/ 

This expression gives the nature of the source function in a 

vibrationally relaxing gas. For X small (high pressure) this 

expression for J(v>z) reduces to the Planck emission B(V,z), 

When an energy source or sink Q(z) outside of the range Av is 


(C.7) 
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present, in view of eqxiation (C«3), the source function iriay be 
expressed as 


J(v,z) = B(v,z) 


1 


Q( z) 


^ ^B(V ,z)n(z)k(V,z)d'i>du) 


(C.S) 


This equation suggests that when the gas is relaxing, it does not 
radiate with Planck intensity in the band, but only with a frac- 
tion of B(v,z), the fraction being given by the term in the 


brackets. 

In order to calculate the temperature structure in atmos- 
pheric regions where radiative transfer of energy is important, 
we assume that radiative equilibrium conditions prevail, so that 
the energy balance described by equation (C.4) holds. Substi- 
tuting from equation (C.B) in (C.4) we have 


Q(z) + j^I(v,e,z)n(z)k(v,z)dvduj 


/■ 


B(V,z)(l-^ 


Q(z) 


9yy'B(v,z)n(z)k('«) ,z)dvdo)^ 


(C.9) 

|n(z)k(v ,z)dvdu; 


or. 


Q(z) + 


JJ I(v ,e,z)n(z)k(v,z)dvdu> 

jj B(v ,z)n(z)k(v ,z)dvdu> 

^ yr/ B(V,z)n(z)k(v,z) 

$ \ ^B(V ,z)n(z)k(v,z)dvdo> 

^B(v,z)n(z)k(v,z)d>>dw is not a function of v or o) but 


(C.IO) 


dudoo 


Since 

only o7z, it may be extracted from the outer integral in the 
second term on the right hand side of (C.IO), and cancelled with 
the numerator. Upon rearrangement of the remaining terms, we 


obtain 
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^^^^'^»®»2)n(z)k(v),z)dydca 
= /b(v f z)n 


(z)k(v ,z)dvdu> 


(C.ll) 


The factor (1 + in the first term on the left hand side of 
(C.ll) may be regarded as a factor amplifying external sources 
and sinks of energy to take into account the fact that the radi- 
ation emanating from a relaxing atmospheric layer is not related 
directly to or characteristic of the kinetic temperature of the 
layer alone. With the source fxmction related to the tempera tiire 
via eq\iation (C.S), it is possible to calculate the radiative 
equilibrium temperature structure of the atmosphere by a simul- 
taneous solution of (C.8) and (C.ll). 

The infrared flux divergence JJ dvdco (equation 
(C.l)) determines the rate at which energy is exchanged between 
the vibrational - rotational energy levels and the radiation 
field. If these energy levels are perturbed, so that a Boltzmann 
distribution is no longer present, then the rate of energy trans- 
fer between vibrational - rotational and translational modes by 
collisions is given by (Curtis and Goody, 1956) 

dE E / E \ , ^ 

»_( 1 - — ) , (C.12) 

dt X \ E / 

where E is the total vibrational energy expressed in quanta per 
unit volume, and E is the value of £ under Boltzmann conditions. 
The "non -Boltzmann” source function J(y,z) will then be related 
to the Planck function by 


£ 

J(v,z) - B(v,z) — 

£ 


(C.13) 
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Equation (C,12) indicates that in high pressure regions, any 
perturbation in the vibrational populations from a Boltzmann 
distribution will rapidly decay. V/hen pressure is low, hov/ever, 
significant departures from a Boltzmann distribution may persist 
indefinitely if, for example, a periodic perturbation such as 
solar heating is operating. But with cessation of a perturbing 
influence, the atmosphere will eventually reach a state of LTE, 
regardless of the collision frequency, as may be inferred also 
from eqxaation (C.ll) if Q(z) is zero. 


APPENDIX D 

TRANSfflSSICN FUNCTIONS OF CARBCN DIOXIDE 
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Line shape and mean line absorption # 

The width and shape of a spectral line are influenced by 
three broadening mechanisms: collisional, Doppler and natural 

broadening. Collisional broadening arises out of molecular 
interactions while Doppler effects are due to molecular ”line-of- 
sight” velocities. Natural broadening results from radiation 
dancing. 

In the case of non-interacting, stationary molecules, 
spectral lines have a finite width due to the uncertainty prin- 
ciple which predicts an xmcertainty AE in an energy level. This 
uncertainty is related to the natural lifetime 9 of the molecule 
in the excited level by the expression (Goody, 1964) 

21T AE 9 = h (D.l) 

where h is Planck’s constant. Hence, all spectral lines invol- 
ving transitions between ground and excited levels will have an 
uncertainty in frequency AV, or natxiral half-width oCj^, (AV/2), 


given by 


Av » AE/h 


or, 



1 

4tT9 


(D.2) 

(D.3) 


The variation of absorption coefficient within a naturally broad- 
ened line may be expressed as (Mitchell and Zemansky, 1934) 

1 


k(V) 


(D.4) 



144 


Thus, natural line broadening is independent of all effects 
external to the molecule. 

Normally, hen*/ ever, natural broadening of spectral lines is 
negligible in comparison v/ith the Doppler effect, resulting from 
the thermal motions of the molecules and therefore related to the 
temperature. The Doppler half-v/idth and line shape are 
readily derived from kinetic theory (Aller, 1953)* The theore- 
tical expression for the Doppler half -width is 

o /2KT\ 1/2 . . 

“d ■ - ( •) . 

c \ m / 


where K is Boltzmann’s constant and m is the mass of the absor- 
bing molecule. We may therefore express the Doppler half -width as 





(D.6) 


where is the value of at reference temperature To* The 
absorption coefficient within a Doppler broadened line is given by 


’wdAT ”'d°) ) ’ 

where S is the integrated line intensity. Doppler broadening 
produces noticeable effects when absorption or emission of a gas 
under low pressures such as those encountered in the higher levels 
of a planetary atmosphere is considered. 

When pressures are high, the line shape is determined by 
collisional broadening. Molecular interactions cause perturba- 
tions of the energy levels of the molecules, and, as a result 
of these perturbations , absorbed or emitted energies are spread 


I 
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over a finite frequency range. In addition, the line center may 
be displaced toward lower frequencies by an amount proportional 
to the total pressure and the line shape may become asymmetrical. 
The classical theory of line shape and absorption coefficient 
within a pressure broadened line was developed by Lorentz, who 


found that SoIt 1 

k(v) = — ^ 5 : 

(V-V„) + a ‘ 


(D.8) 


where is the half -width of the line due to collisions. The 
theoretical evaluation of the Lorentz half -width is difficult, 
and experimental values are normally used. From kinetic theory, 
the collision frequency depends upon the ratio p//T, so that 
which is proportional to the collision frequency, is given by 

‘^L ' , (3.9) 

where WLoI® measured Lorentz half-v/idth at reference pressure 
p® and temperature To, and p^= is an effective broadening pressure. 
For collisional broadening of CO 2 lines, p*!* may be expressed as 

p’^ = p + .3 pc02 » (D.IO) 

where PCO 2 partial pressure of C02» In a pure CO 2 atmos- 

phere, a correction of about 30^ to the total pressure is re- 
quired to properly take into account the more effective self- 
broadening of CO 2 in computing the Lorentz half -width. 

In addition to the Lorentz theory of collisional line shapes 
a number of alternate theories and modifications have been pro- 
posed (cf. Lindholm, 1945; Benedict, ^ al« , 1956). There is 
nevertheless strong support for use of the basic Lorentz line 


profile (Curtis and Goody, 1956; Young, 1964), and it appears 



that, due to the other inaccuracies involved in a study of this 
type, it would serve no purpose to consider non-Lorentsian line 
shapes for collisional broadening in lov;er layers of an atmosphere. 

V/hen all three broadening mechanisms act togetlier, the line 
absorption coefficient takes the form (Mitchell and Zemansky, 1934) 



dll 

(V -Vo + UVo/c)^ + 


(D.ll) 


or, in view of (D.5), 
k(V) = S 



ys-J) cjtp/ _Vc U ^ 2 cx 


M- 


^ (V-^4+ UVo/c)^'^ +06^^ 


(D.12) 


where U is the mean molecular velocity. Taking the substitutions 

/ mlJ^ \ 1/2 ,n 

y \ 2XT / » (D»13) 

/. .A 

. (D.14) 


V = 




k(0) = 


and 


ot-Q ’ 


CX=(X,/o(„ , 


(D.15) 

(D.16) 


equation (D.12) may be written (Plass and Fivel, 1953) 


-VOO 


= — / — ^ dy . 

k(0) 1V/ <x2 + (v-y'; 


(D.17) 


Under low pressures, (<^<.5)i Harris (194^)) has shov;n that the 
right-hand side of (D.17) can be expressed as a Taylor series, 
H(cx ,v) , so that 
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Band nodels . 

:^quations (D»20) and (D*21) represent the r.ean fractional 
absorption associated v;ith a single spectral line. The calcu- 
lation of the absorption or transmission over a finite spectral 
interval is extremely difficult due to the larger number and 
complex arrangement of the rotational lines in a vibration - 
rotation band. Various models have been proposed to represent 
the mean transmission function of the lines in a real band. The 
"regular" or Elsasser model (Elsasser, 193^, 1942) is based upon 
the assumption that the lines all have the sane intensity and 
half-v/ifth and are uniformly spaced. The "statistical" nodel 
(Goody, 1952 ) assumes a random spacing of the line centers v/itli 
the intensities represented by a probability function. The 
"random Elsasser" model of Plass (195S) is somexvhat less restrict- 
ive, with the spectrum represented by a random distribution of 
several Elsasser bands, each with its characteristic line inten- 
sity and line spacing. Recently, VJyatt, ^ e_l, (1962) have 
introduced the "quasi-random" model, in which a spectral interval 
of interest is divided into a number of much smaller intervals, 
each of which contains a random distribution of lines. The lines 
in each sub-interval are grouped according to intensities, and 
the contribution to the absorption in the interval from the wings 
of lines outside of the interval is taken into account. This 
model is probably the most realistic representation of a spectral 
band introduced to date. 
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v;here the coefficients Hj^(v) are tabulated functions. The 
absorption coefficient v/ithin the line is thus given by 


k(v) =^H(o<,v) . (D.19) 

From this expression for the absorption coefficient, the mean 
fractional absorption under low pressures for a spectral line, 
over an interval of width d equal to the mean line spacing, can 
be expressed as 

0( /'*d/2«D 

X » - exp(-k(v)u)) dv (D.20) 

J -d/20Cu ^ ^ 

v/here u is the optical path. "K is evaluated numerically by 
substituting the tabulated functions Hj^(v) of Harris (194^) in 
(D.19) for k(v) . 

Van der Held (1931) has shown that when pressures are high 
(<^>1), the mean fractional absorption for a spectral line of 
Lorentz shape is closely approximated by the Ladenburg-Reiche 
( 1913 ) formula 

r “ 1 - exp(^^xe*’^(Io(x) + l3 ^(x) )^ , (D.21) 

where Iq(x) and I]_(x) are the Bessel functions of order zero 
and one and of imaginary argument, and v/here 


/3 = 



(D.22) 


and 


X 


Su 




(D.23) 
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For calculations of temperature from energy balance consid- 
erations, it can be seen that, if the dominant terms in eouation 
(9) are those representing the energy absorbed and emitted per 
unit volume in the 15 CC 2 band, the effects of the local trans- 
mission function in this band on the temperature will be minimal, 
since this function has equal v;eight in both terms. Due to this 
compensatory effect, the choice of a particular band model does 
not critically affect the determination of atmospheric 
temperature* However, care must certainly be exercised in the 
choice of a band model when time dependent calculations involv- 
ing the evaluation of beating rates are performed. For present 
purposes, the transmission or absorption of a gas can be 
adequately described by any one of several band models. For a 
linear, "well behaved" molecule such as C 02 > line intensities and 
spacing are fairly uniform, and the Slsasser model closely 
approximates the true absorption (Kov/ard, ^ al, 1955). 

In the present study, a modified version of tr^e Flsasser 
model is used. The band of interest is divided into several 
frequency intervals. The Lorentz half-width and mean line spac- 
ing are assumed to be the same throughout the band, while the 
Doppler half-v/idth and mean line intensity are varied from 
interval to interval within the band. For such a model, the 
mean fractional absorption for each spectral interval becomes 
identical v/ith the mean fractional absorption for an individual 
spectral line given by equation (D.20) or (D.21), 





In a planetary atmosphere, pressure varies by many oraers 
magnitude in the vertical. Variations of temperature are ordin- 
arily much smaller, but these may also be significant. In using 
the above band model to evaluate transmission functions over 
long atmospheric paths, the pressure and temperature variations 
must be taken into account. Since the band parameters S, and 
oi. Q vary with pressxire and temperature, equations (D.7) and (D.S) 
cannot be used to obtain the absorption coefficients for an atmos- 
pheric path over which p and T variations are significant. 

Curtis ( 1952 ) and Godson (1953) derived an approximate method 
for dealing with variations of pressure along an absorbing path. 
Their method consists essentially in choosing a weighted value 
for pressure along the path which, when used to evaluate the 
fractional absorption from (D.20) or (D.21), results in the 
same value of A as would be obtained by integration of the 
absorption coefficient over the path. For isothermal paths 
along which the CO 2 mixing ratio is constant, the weighted pres- 
sure given by the Curtis-Godson approximation is simply the 
average of tne pressures at the two end points. The mean tem- 
perature for the path may also be taken to be the average of the 
two end points. Unless large temperature maxima or minima occur 
in the atmosphere, such a mean value should be fairly character- 
istic of the temperature along the path. 

The Curtis-Godson approximation has been found satisfactory 
for transmissivity and flux calculations (Godson, 1953 J Kaplan, 
1959 ) and also for heating rate calculations (Walshaw and 
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Rodgers, 1963 )♦ It was used in the present study to evaluate 
the parameter o( for any atmospheric path required. 

Using the Curtis-Godson approximation together with the band 
model described above, the transmission functions were evaluated 
from equations (D.20) and (D.21). In the lower levels of a 
planetary atmosphere, where the pressure is relatively high, an 
appropriate form for A is the Ladenburg-Reiche expression (D.21). 
In the higher, more rarified atmospheric regions, equation (D.20), 
based on Harris (1946), properly represents the true absorption, 
while at intermediate levels, a suitable interpolation scheme 
may be employed. 


Band parameters . 

Stull, ^ ( 1963 ) have made a detailed theoretical inves- 

tigation of the transmission of CO 2 in the spectral region 
between 1 ytf and 7 . 0 . They have obtained the absorptive prop- 
erties of COg for effective broadening pressures ranging from 
1 atm to 10 “ atm and path lengths from 10^ cm-atm to 0.2 cra-atm 

for three temperatures: 200*K, 250*K and 300*K. All contributing 

-6 

spectral lines whose relative intensity is greater than 10 
times that of the strongest line in any particular band were 
included, and the various isotopes of C and 0 were taken into 
account. Transmission functions were calculated using the 


"quasi-random” band model and assximing a Benedict line shape for 
refined spectral intervals, as small as 2.5 cm”^. In Figure D.I, 
the absorption spectrum of COg in the lyt{ - 6 ^ near infrared 
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region obtained by Stull, et (1963) is shovm for optical 
4 2 

paths of 10 cra-atm and 10 cm-atm. Comparison of the theoreti- 
cal transmissions of Stull, at al, (1963) with the measurements 
of Burch, at al. (1962) in tha strong CO 2 bands at 2.7ytf , 4.3^ 
and 15 ^ indicatas axcellent agraement. 

Estimatas of tha band parameters required for the evaluation 
of the near and far infrared transmission functions of CO^ v/ere 
obtained from the transmission tables of Stull, ^ (1963) in 

the following manner. 

For constant values of ^ (D.22) and constant values of 
(D.23), a log-log plot was made of the Ladenburg-Reiche function 

f^(^,x) “ 1 - exp^^x e“^ (Iq(x) I^(x))^ (D.24) 

against the parameters ^^x. Also plotted against ^x was the 
function 

- 1 - . (D. 25 ) 

The function f^(^,x) is the exact expression for the mean absorp- 
tion in a spectral interval given by the statistical model, and 
the function f^(^,x) is the absorption given by the "strong line" 
approximation to the statistical model, applicable for large 
values of/3x and small values of^ (Plass, i 960 ). The network 
for f^(^,x) shown in Figure D.II was thus obtained, v/ith the 
uppermost line in the network corresponding to The 

functions f^(^,x) and 

=» 1 - exp(-^x) 


(D.26) 




Kure U.II. (/3 ,x) vs /6'-x for constant values of /3 agd 
” constant values of /5x, Also f2(./6»^) vs x 
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were also plotted in log-log fashion against the parameter 

at constant values of ^(Figure D.III)* f^(^,x) represents 
the mean absorption in a spectral interval as given by the ’’weak 
line" approximation, valid when j3 is large and x is small (Plass, 
I960) . 

From the tables presented by Stull, et (1963), the 
average fractional absorption A for a frequency interval of inter- 
est was extracted for various path lengths u and pressures p* at 
a temperature of 200*K. The mean absorption at constant u and 
constant p was plotted in log-log fashion a^ a function of the 
parameter up , and the mean absorption at constant p was plotted 
in the same way as a function of u. As an example, plots of A 
at constant u and p* against up and of A at constant p* against 
u for the far infrared interval between 710 cm"^ and 730 cra“^ 
are presented in Figures D.IV and D.V. 

The relationships between ^ x and up and ^x and u are evi- 
dent from equations (D,22) and (D«23)» By superimposing the 
patterns of mean fractional absorption in a spectral interval 
from Stull, ^ (1963) (such as Figures D.IV and D.V) on the 

general patterns for A (Figures D.II and D.III), it is possible 

2 

to obtain a correspondence between ^ x and up ,^x and u and ^ 

and p , In spectral intervals where the absorption is strongest, 

only a "strong line” fit could be obtained, giving the relation- 
2 ^ 

ship between^ X and up • In weaker intervals, however, it was 
generally possible to obtain the correspondence between^^x and u 
and ^ and p also. Thus, in some spectral intervals, a consist- 
ency check on the estimated values of the band parameters was 
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provided. The ratios 


2irdL,So 

' ■ B “■ ■■ 


up 

d^ 

/3x 

A 

u 

d 

A « 


P* 

d 


where the subscript (o) refers 


, (D.27) 

, (D.2S) 

, (D.29) 

to 200® K, may be manipulated to 


yield values for the parameters 



(D.30) 



{D.31) 


(Thus X “ C u/p* and /6“ C p . ) The preliminary estimates of 
ill 2 

and obtained by curve fitting were corrected by trial and 
error to find the best possible agreement with Stull, ^ al. 
(1963). Final estimates of and at a temperature of 
273. 2®K for various spectral intervals in the near and far infra 
red are listed in Table D.I. 

Yamamoto and Sasamori (1956) have used a mean line spacing 
d of 1.55 cm“^ in their calculation of the absorption in the 
15 /f CO 2 band. The mean half -width for pressure broadening in 
this band C<t has been determined by Kaplan and Eggers (1956) 
usins a curve of growth method as O.O 64 cm at 1 atm and 296® K. 

Correcting this estimate of 0(^ to 273. 2*K, a value for of 
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TABLE D.I 


BAND PARAMETERS C AND PCJR NEAR AND FAR INFRARED CO BANDS 

12 2 


Band Interval (cm“^) 


(atm cm”^) 


C (atm*”^) 
2 


1.3« 7775-7725 

' 7725-7675 

7675-7625 
7625-7575 
7575-7525 
7525-7475 
7475-7425 
7425-7375 

1.4/1 7025-6975 

' 6975-6925 

6925-6S75 

1.6/1 6575-6525 

' 6525-6475 

6475-6425 
6425-6375 
6375-6325 
6325-6275 
6275-6225 
6225-6175 
6175-6125 
6125-6075 
6075-6025 

2 . 0/1 5200-5150 

' 5150-5100 

5100-5050 
5050-5000 
5000-4950 
4950-4900 
4900-4650 
4650-4600 
4600-4750 
4750-4700 

2 .7/1 3650-3600 

' 3600-3750 

3750-3700 

3700-3650 

3650-3600 


5.0 X 

10“^ 

.30 

1.2 X 

10“5- 

.30 

1.5 X 

10-5 

10-f 

.30 

2.0 X 

.30 

1.3 X 

10-4 

.30 

2.0 X 

10 “i*- 

.30 

7.2 X 

10"? 

10"5 

.30 

1.6 X 

.30 

9.2 X 

000 

1 1 1 

.30 

2,0 X 

.30 

1.2 X 

.30 

2.9 X 

10"t 

.30 

9.0 X 

.30 

2.9 X 
6.4 X 

10 ? 
10“? 

.30 

.30 

7.0 X 

10"^ 

.30 

2.5 X 

10"5 

.30 

3.5 X 

10”^ 

.30 

2.3 X 

10’^ 

.30 

1.0 X 
5.5 X 

10:^ 

10 5 

.30 

.30 

1.9 X 

.30 

3.9 X 

10"^ 

.33 

3.5 X 

10 ~l 

.33 

3.0 X 

10”^ 

10"3 

.33 

7.5 X 

.33 

5.3 X 

10-1 

.33 

2,0 X 

lOll 

.33 

1.0 X 

6.0 X 

10 i 
10", 
10“ f 

.33 

.33 

1.3 X 

.33 

3.0 X 

10-^ 

.33 

1.0 X 

lo"!^- 

.36 

2.0 X 

10-J 

.36 

2.7 X 

10 1 

.36 

5.0 X 

10 ? 

.36 

1.7 X 

10 1 

.36 


(Continued) 
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TABLE D,I (Continued) 


BAND PARAI*ETERS C, AND C 

1 2 


FOR NEAR AND FAR INFRARED CO 

2 


BANDS 


Band Interval (cm**^) 


C^ (atm cm”^) 


C^ (atm“^) 





3600-3550 

6»0 

X 

10 ? 

.36 

3550-3500 

2.1 

X 

lO'i- 

.36 

3500-3450 

5.0 

X 

10-3 

.36 

3450-3425 

2,0 

X 

10-5 

.36 

2470-2450 

2.2 

X 

10“5 

.66 

2450-2430 

2.0 

X 

10-5 

.66 

2430-2410 

5.0 

X 

10-^ 

.66 

2410-2390 

3.2 

X 

10-3 

.66 

2390-2370 

1.6 

X 

10 ? 

.66 

2370-2350 

7.0 

X 

10 t 

.66 

2350-2300 

2,0 

X 

10 t 

.66 

2300-2250 

7.0 

X 

10-1 

.66 

2250-2200 

6.0 

X 

10-3 

. 66 

2200-2150 

6.6 

X 

10-5 

.66 

2150-2100 

1.0 

X 

10-3 

.35 

2100-2050 

2.0 

X 

io~l 

.35 

2050-2000 

1.0 

X 

10“3 

.35 

2000-1950 

5.0 

X 

10-5 

.35 

1950-1900 

2.0 

X 

10-3 

.35 

750-730 

1.1 

X 

10-^ 

.27 

730-710 

5.6 

X 


.27 

710-690 

2.4 

X 

10 0 

.27 

7OO-6SO 

S.5 

X 

10 ? 

.27 

690-670 

3.6 

X 

10 } 

.27 

6S0-660 

9.9 

X 

10 t 

.27 

670-650 

9.9 

X 

10 1 

.27 

660-640 

2.4 

X 

10 n 

.27 

650-630 

6.S 

X 

10 0 

.27 

640-620 

2.0 

X 

10 ? 

.27 

630-610 

6.3 

X 

10-1 

.27 

620-600 

4.4 

X 

10-1 

.27 

SOO-750 

7.6 

X 

10-? 

.27 

750-700 

4.1 

X 

10-1 

.27 

700-650 

3.6 

X 

10 i 

.27 

650-600 

2.7 

X 

10 0 

.27 

600-550 

1.3 

X 

10-2 

.27 
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0,27 at this temperature is obtained. Also, the mean line inten- 
sity for a chosen interval within the band can be estimated, 

in view of this knowledge of d, f from equation (D«30), 

Although they are not directly involved in the calculation of the 
transmission functions, estimates of S (at 200* K) were obtained 
and are presented in Table D,II. 

No estimates of 0(^ or d are available for the - 6^ 
near infrared CO^ bands. However, as can be seen from Table D.I, 

the estimates of C for the near infrared bands are quite close 

2 

to that for the 15 fi region, the largest discrepancy being found 

in the 4.3yt/ band, where is larger than the 15^ value by a 

factor slightly greater than three. This implies that both (X^ 

and d are of approximately the same size in both the 1 ji - 6^ 

region and 15yH region. Then, if the Kaplan and Eggers (1956) 

value of 0{. is adopted in the near infrared bands also, an 
L 

estimate of d in each near infrared interval can be derived from 
(D.31) and S can be obtained from (D.30). The values of d and 
S (at 200*K) resulting for spectral intervals in the 1 j4 - 6 
region are given in Table D,II, 

The temperature dependence of the mean line intensity for 
each spectral interval in the 15 ji CO^ band can be taken into 
account by an exponential term representing the increase of the 
population of the lov/er vibrational energy state with temperature 
(Burch, et al, . 1961). This temperature dependence in the 1 - 
6yt{ bands is negligible, as is a small effect due to the temper- 
ature dependence of the rotational levels associated with each 
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TABLE D.II 


BAND PARAl'^IETERS S, d A>iD P FOR INFRARED CO^ BANDS 


Band Interval (cm“^) 


1.3 M 7775-7725 

' 7725-7675 

7675-7625 

7625-7575 

7575-7525 

7525-7475 

7475-7425 

7425-7375 

1.4w 7025-6975 

6975-6925 

6925-6675 

1.6 m 6575-6525 

6525-6475 
6475-6425 
6425-6375 
6375-6325 
6325-6275 
6275-6225 
6225-6175 
6175-6125 
6125-6075 
6075-6025 

2 . 0/4 5200-5150 

' 5150-5100 

5100-5050 
5050-5000 
5000-4950 
4950-4900 
4900-4650 
4650-4600 
4600-4750 

4750-4700 


-1 -2 -1 

S (atm cm ) d (cm ) 


2,1 

5.0 

X 

X 

10“i 

l°l6 

1.39 

1.39 

6.3 

6.4 

X 

X 

107 

10-^ 

1.39 

1.39 

5.3 

X 


1.39 

6,4 

X 

10 A 

1.39 

3.0 

X 

10“^ 

1.39 

7.5 

X 

10“° 

1.39 

3.6 

X 

10"^ 

1.39 

6.4 

X 

10“^ 

1.39 

5.0 

X 

10“^ 

1.39 

1.2 

X 

10“^ 

1.39 

3.7 

X 


1.39 

1.2 

X 

10 1 

1.39 

2.7 

X 

10“° 

1.39 

2.9 

X 

10“J 

1.39 

1.0 

1.5 

X 

X 

10“^ 

1.39 

1.39 

9.6 

X 

lOll 

1.39 

4.2 

X 

10 ° 

1.39 

2.3 

X 

loil 

1.39 

7.9 

X 

10 ° 

1.39 

1.6 

X 

10"^ 

1.27 

1.5 

X 

10“2 

1.27 

1.3 

X 

10“? 

1.27 

3.1 

X 

10"? 

1.27 

2.2 

X 

10“1 

1.27 

6,4 

X 

10“3 

1.27 

4.2 

X 

10“2 

1.27 

2.5 

X 

10“f 

1.27 

5.4 

X 

10“^^ 

1.27 

1.3 

X 

10"^ 

1.27 


rCK ) 


(Continued) 
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TABLE D.II (Continued) 


BAND PAR;U-IETERS S, d AND P FOR INFRARED CO2 BANDS 


Band Interval (cm“^) S (atm“^cm"^) d (cm“^) 


P (*K) 


2 



4 . 3 ^ 



3650-3600 

4.2 

X 

10-5 

lo-;*- 

3600-3750 

6.4 

X 

3750-3700 

1.1 

X 

0 

10 0 

3700-3650 

2.1 

X 

3650-3600 

7.1 

X 

10 g 

3600-3550 

2.5 

X 

10 9 

3550-3500 

8.8 

X 


3500-3450 

2.1 

X 


3450-3425 

8.4 

X 

10 ^ 

2470-2450 

9.2 

X 


2450-2430 

6.4 

X 

10 ° 

2430-2410 

2.1 

X 

10 “^ 

2410-2390 

1.3 

X 

10-f 

10-1 

2390-2370 

6.7 

X 

2370-2350 

2.9 

X 

0 

2350-2300 

6.4 

X 

10 ^ 

2300-2250 

2.9 

X 

10 

2250-2200 

2.5 

X 

10 ^ 
10 ^ 

2200-2150 

2.7 

X 

2150-2100 

4.2 

X 

2100-2050 

6.4 

X 

10-3 

2050-2000 

4.2 

X 

10“2 

2000-1950 

2.1 

X 

10-5 

10“4 

1950-1900 

6.4 

X 

750-730 

4.6 

X 

10“2 

730-710 

2.3 

X 

10-1 

710-690 

1.0 

X 

10 0 

700-680 

3.6 

X 

10 0 

690-670 

1.6 

X 

10 1 

68O-66O 

4.1 

X 

10 1 

670-650 

4.1 

X 

10 1 

66O-64O 

1.0 

X 

10 0 
10 0 

650-630 

3.7 

X 

640-620 

6.4 

X 

lOIi 

630-610 

3.5 

X 

10*1 

620-600 

1,6 X lO"-^ 

( Continued ) 


1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

.465 

.465 

.465 

.465 

.465 

.465 

.465 

.465 

.465 

.465 

1.19 

1.19 

1.19 

1.19 

1.19 

1.55 

1.55 

1.55 

1.55 

1.55 

1.55 

1.55 

1.55 

1.55 

1.55 

1.55 

1.55 
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TABLE D.II (Continued) 



BAI\iD PARiU'STERS S, d 

AND r FOR INFRARED 

CO 2 BANDS 


Band 

Int erval s ( cm"^ ) 

S (atm”^cra”^) d 

( cm"^ ) 

r (‘'K) 

15/< 

BOO-750 

3.2 X 10"? 

1.55 

2.0 X 10? 


750-700 

1.7 X 10"f- 

1.55 

1.3 X 103 


700-650 

1.6 X 10 1 

1.55 

4.0 X 10^ 


650-600 

1.1 X 10 ? 

1.55 

1.0 X 10| 


600-550 

5.4 X 10-3 

1.55 

2.0 X 10^ 
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vibrational state in the 15 ji region, v/hich alters the shape of 
individual bands (Hanel and Bartko, I 964 ). The mean line inten- 
sity in a particular spectral interval of the 15 yij band may 
therefore be expressed as 

3(T) = S, exp(r (t.“^ - , (D. 32 ) 

where P is the coefficient of temperature dependence and S® is 
the mean line intensity at temperature T« » Thus, from a know- 
ledge of S at 200*K, S at any other temperature may be computed, 
if r is knovm. 

The parameter P in a sub-interval of the 15 ji band v;as 
obtained from equation (D,32) using the transmission functions 
of Stull, ^ a^. ( 1963 ) at three different temperatures. 

Obtaining S at 200* K, 250*K and 300* K by the procedure outlined 
above, values of P , representing an empirical fit to equation 
(D. 32 ) were derived for selected spectral intervals. The esti- 
mates obtained for the coefficient of temperature dependence 
for various 15/( band intervals are given in Table D.II. 

Combined with estimates of the other band parameters C^, 

Cp, d, (X , (X and 3*, the determination of P enables us to 

li 0 

eva.luate the atmospheric transmission functions required for 
solution of the radiative transfer equation from expressions 
(D.20) and (D.21). In the weaker bands of the 1^ - jA region, 
these estimates of the various band parameters lead to precise 
agreement with the transmission functions of Stull, et (1963). 
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Even in the strongest spectral regions near 4.3yi< and 15^ , a 

comparison of the transmission functions obtained here with 
those of Stull, indicates close agreement (Figure D.VI). 

Numerical evaluation of COo transmission functions . 

The determination of atmospheric temperatur’es by a sceady- 
state, iterative solution of the radiative transfer equation 
requires the evaluation of tens of thousands of transmission 
functions. Direct evaluation of a making use of equations (D.20) 
and (D.21) is quite involved, and innumerable repetitions of this 
procedure are clearly not practical, even making use of the most 
modern high speed electronic computers. Therefore, a more eco- 
noieical method was developed, utilizing a table-searching and 
interpolation scheme to obtain the necessary transmission 
functions. 

Using the value of at the center of the 15^ CO^ band 
( 667.40 cm and at a temperature of 150*K, along with the esti- 
mate of d given by Yamamoto and Sasamori (1950), a table of the 
parameter (X d/ Oi ) was constructed, with coordinates oc and 

2^ OCX, from equations (D.20) (for^^<l) and (D.21) (forex >1). 

—6 2 

Enti’ies were made for values of ranging from 10“ to 10 in 

order-of-magnitude steps, and for values of T'iVax ranging from 

—1 9 

10 to 10 in steps of tenths of orders-of-r.agr.itU'V-, Thus, 
this table has dimensions 9x91, v/ith a total of Sl9 entries. A 


plot of these entries is shovm in Figure D.VII, 
tional line has been added corresponding to o>c 








Figure D.VII represents X as a function of ^ arid 2<ra;: at 
667.40 cn~^ and 150”K; a slightly modified version v.’ould be 
representative of absorption at another v/avonuniber and/or temper- 
ature, v;ith differences occuring tovjard larger values of 
(saturation) for a given oL , ‘n upward or dovmward sViift in the 
position of the saturation level (uppermost curve) and in the 
location of lines of constant ^ near this level would result 
from consideration of frequencies and temperatures other than 
those for vniich the figure v/os constructed. All except the 
upperitost portion of Figure D.VII 'v/ould, hov;ever, remain 
unchanged, unless drastically different values of fx’equency and 
temperature v;ei'e used. Therefore, a suitable corx’oction scneme 
was .adooted to adjust A near the saturation leveu ior frequency 
and tenperatxire , so that a single table could oe used to ootain 
the absorotion for all of the frequency intervals in the u5 ji 
band at any desired temperature. 

For a given atmospheric path the parameters 

and 2<botx (=Su/ ) were evaluated for a particular frequency 

D 

interval using the weighting procedure described previously. 

Then, by logarithKiic interpolation on the cocruinate and a 
linear interpolation on the 2'tbodx coordinate, a value of 
(X d/ a )'*' was obtained from the table, where the asterisk corre- 
spends to 667.40 cm“^ and 150*K. (This interpolation method was 
found to reproduce quite accurately the values of (A d/ a^) 
computed directly for intermediate values of (X and 2^<xx from 
equations (D.20) and (D.21).) The value of (A d/Of,- ) ’ extracted 
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in this way from the table vra.s then corrected for the variation 
of (X^ with frequency and/or temperature to obtain 3T, making 
use of the expression 



(D.33) 

(This formulation leads to values of T in excellent agreement 

with those calculated directly from (D.20) and (D.21) taking 

the frequency and temperature dependence of Oi into account*) 

D 

A scheme was also devised for making use of the table of 
(7T d/o^^) to obtain 3T whenever the values of oi and 2'lvax fall 
outside of the covered range. In the Martian atmosphere, with 
a surface pressure near 10 mb, the value of for any path can 
never be much greater than 10^, and must certainly be less than 

2 . .4 

10 ; it may, on the other hand, become smaller than lo” • V/hen- 
ever a is this small, however, the value of 2'frax must be rela- 
tively small also. Thus, when oi lO"^, (X d/<X^)*^ for a 
given 21V<xx will be independent of ot ^ and equal to the value 
of (X d/oljj)’'' at = 10"^, as indicated by Figure D.VII. V/hen 
^ < 10“^, all required values of (X d/oij^)^ fall on the 

lowermost curve in the figure, regardless of the particular 
value of Of , and X is determined by 2*<Vocx alone. For some path- 
lengths in the atmosphere of Mars, the parameter 2^otx lies 
outside the range of the table also. However, when 2lV<xx is 
smaller than 10~‘*‘, (A d/ 0i_) is equal to 2'(Vax for all <x (Fig- 

^ 9 

ure D.VII). In addition, whenever 2trax exceeds 10 , ot too 
must be relatively large, so that saturation has been reached 
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and { A d/ 0^ ) assumes its maximum value, V/ith these considera- 

^ __ Jjc 

tions, the tabulated values of (A d/ Oi^) could be used to deter- 
mine A for any required path, throughout the entire range of 
temperature, pressure and optical thickness encountered in the 
Martian atmosphere, in any chosen spectral interval in the 15^ 
band. 


A similar procedxire v/as used to evaluate A for spectral 
intervals in the 1^ - 6 ^ CO 2 bands. Due to its frequency depend 
ence, the Doppler half-v/idth in these bands is larger by nearly 
an order of magnitude than 0^^ in the 15y^f band. Figure D.VII, 
therefore, would have to he considerably modified to properly 
represent A in these bands, with saturation occuring at a lower 
level, and correction of the tabular values of (A d/o^j^)’‘‘, 
using equation (D.33) would be of doubtful accuracy for some 
atmospheric paths. Therefore, to avoid the lengthy evaluation 
of A via (D.20) or (D.21), another table of the same dimensions 
as that described previously v/as constructed, using the value 
of near the center of the 1,6^ region (6325 cm ) at a 
temperature of 150*K, and the estimated value of d in the 1.6^ 
band given in Table D.II. Operating on this second table in 
exactly the same manner as described above, and using tne same 
correction scheme (equation (D,33))» fractional absorption in the 
desired intervals in the near infrared was obtained. 

The near infrared transmission function for various paths 
obtained from the table were compared with those calculated 
directly from (D,20) and (D,21) and the agreement was founa to 
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be quxte satisxactory. This method of obtaining transmission 
functions in the - 6^ CO^ bands was used in determining 
the transmission of near infrared solar energy in the Martian 
atmosphere. 
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